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I. INTRODUCTION 



Black holes with a wide range of masses could have formed in the early Universe as a result of the great compression 
associated with the big bang [lJ-Q . A comparison of the cosmological density at a time t after the big bang with the 
density associated with a black hole of mass M suggests that such "primordial" black holes (PBHs) would have a 
mass of order 

This roughly corresponds to the particle horizon mass in a noninflationary model or the Hubble mass otherwise. PBHs 
could thus span an enormous mass range: those formed at the Planck time (lCU 43 s) would have the Planck mass 
(10 _5 g), whereas those formed at Is would be as large as 10 5 Mq , comparable to the mass of the holes thought to 
reside in galactic nuclei. They could be even larger than this in some circumstances. By contrast, black holes forming 
at the present epoch could never be smaller than about 1 Mq . 

The high density of the early Universe is a necessary but not sufficient condition for PBH formation. One possibility 
is that there were large primordial inhomogeneities, so that overdense regions could stop expanding and recollapse [J- 
Hj]. The quantum fluctuations arising in various inflationary scenarios are of particular interest in this context, as has 
been discussed in numerous papers. In some of these scenarios the fluctuations generated by inflation are "blue" (i.e. 
decrease with increasing scale) and this means that the PBHs form shortly after reheating |7Hl(j] . Others involve some 
form of "designer" inflation, in which the power spectrum of the fluctuations — and hence PBH production — peaks on 
some scale |lll425j |. In other scenarios, the fluctuations have a "running index," so that the amplitude increases on 
smaller scales but not according to a simple power law (26l - [33j . Finally, PBH formation may occur due to some sort of 
parametric resonance effect before reheating [34H40l ]. In this case, the fluctuations tend to peak on a scale associated 
with reheating. This is usually very small but several scenarios involve a secondary inflationary phase which boosts 
this scale into the macroscopic domain. 

Whatever the source of the inhomogeneities, PBH formation would be enhanced if there was a sudden reduction in 
the pressure — for example, at the quark-hadron era (4ll443T ] — and especially likely if the early Universe went through 
a dustlike phase at early times as a result of either being dominated by nonrelativistic particles for a period p4l446| or 
undergoing slow reheating after inflation [H, H3| ■ Another possibility is that PBHs might have formed spontaneously 
at some sort of phase transition even if there were no prior inhomo gen eities — for example, from bubble collisions 
[481453 ] or from the collapse of cosmic strings (551463} or necklaces [U, |65| or domain walls [6q-l7l| . 

These PBH formation scenarios are reviewed in Refs. [72|, [73|. Although we do not discuss them in detail here, it 
should be stressed that in most of them the PBH mass spectrum is narrow and centered around the mass given by 
Eq. (jl.II) with t corresponding to the reheating epoch or the time at which the characteristic PBH scale reenters the 
horizon in the inflationary context or to the time of the relevant cosmological phase transition otherwise. However, 
PBHs may be smaller than the horizon size at formation in some circumstances. For example, PBH formation is an 
interesting application of "critical phenomena" (73478} and this suggests that their spectrum could be more extended 
and go well below the horizon mass (79j48l|. This would also apply for PBHs formed during a dustlike phase [8^ |. 
Note that a PBH could not be much larger than the value given by Eq. (jl.II) at formation because it would then be 



a separate closed Universe rather part of our Universe [83j. However, it could still grow subsequently as a result of 



accretion, so the final PBH mass could well be larger than the horizon mass at formation. 

The realization that PBHs might be small prompted Hawking to study their quantum properties. This led to his 
famous discovery [H], HH that black holes radiate thermally with a temperature 
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so they evaporate completely on a time scale 
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T( M)~r_^ 10 ^__j yr . (1 . 3) 

Only PBHs smaller than M* ~ 10 15 g would have evaporated by the present epoch, so Eq. (jl.ip implies that this effect 
could be important only for ones which formed before 10~ 23 s. Since PBHs with a mass of around 10 15 g would be 
producing photons with energy of order 100 MeV at the present epoch, the observational limit on the 7-ray background 
intensity at 100 MeV immediately implied that their density could not exceed about 10 -8 times the critical density 
[86j . This suggested that there was little chance of detecting their final explosive phase at the present epoch, at least 
in the standard model of particle physics 87]. It also meant that PBHs with an extended mass function could provide 
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the dark matter only if the fraction of their mass around 10 15 g were tiny. Nevertheless, it was soon realized that 
the 7-ray background limit does not preclude PBHs having important cosmological effects [88] . These are of different 
types, depending on the PBH mass range. 

PBHs with M > 10 15 g. These would still survive today and might be detectable by their gravitational effects. 
Indeed such PBHs would be obvious dark matter candidates. Since they formed at a time when the Universe was 
radiation dominated, they should be classified as nonbaryonic and so could avoid the constraints on the baryonic 
density associated with cosmological nucleosynthesis. They would also be dynamically cold at the present epoch and so 
would be classified as cold dark matter (CDM). In many respects, they would be like (nonbaryonic) weakly interacting 
massive particles (WIMPs) but they would be much more massive and so could also have the sort of dynamical, lensing, 
and gravitational- wave signatures associated with (baryonic) massive compact halo objects (MACHOs). At one stage 
there seemed to evidence for MACHOs with M ~ 0.5 M from microlensin g ob servations [89| and PBHs formed at 
the quark-hadron phase transition seemed one possible explanation for this [41]. The data now seems less clear but 
there are no constraints exclu ding PBHs in the sublunar range 10 20 g < M < 10 26 g [l9l - [2~l| or intermediate mass 
range 10 2 M Q < M < 1O 4 M [lS|4lj9ifro m having an appreciable density. Large PBHs might also influence the 
development of large-scale structure 9T95J, seed the supe rmassive black holes thought to reside in galactic nuclei 



[9^499j ] . generate background gravit ation al waves [lO0l4lO4 |. or produce x rays through accretion and thereby affect 

the thermal history of the Universe [l05| . 

PBHs with M ~ 10 15 g. As already noted, these would be evaporating today and, since they are dynamically cold, 

one would expect some of them to have clustered within the Galactic halo. Besides contributing to the cosmological 

7-ray bac kground, an effect we reassess in this paper, such P BHs could contribute to the Galactic 7-ray background 
106 . Il07 1 and the antiprotons or positrons in cosmic rays (88l . Il08l . 1 109j ] . They might also generate gamm a-ray bursts 
TToj , radio bursts [m| and the annihilation-line radiation coming from center of the Galaxy [l!2l 113 . The energy 



distribution of the particles emitted could also give significant information about the high-energy physics involved in 
the final explosive phase of black hole evaporation |l!4| . 

PBHs with M < 10 15 g. These would have completely evaporated by now but many processes in the early Universe 
could have been modified by them. F or ex ample, PBH evaporations occurring in the fir st second of the big bang could 
gene r ate t he entropy of t he U niverse [l!5j | , change the details of ba r yoge nesis 1 16Ml 20j | , provi de a s ourc e of neutrinos 



gene i 

[ll3l Il21j or gravitinos [1 22j j or o ther hypothetical particles |123l . |124{ , swallow monopoles [125l . 1 1261 ] , and remove 
domain walls by puncturi ng them [l27l| . If the evaporations left stable Planck-mass relics, these might also contribute 
to the dark matter [jj, l2llll28i - ll36j . Most strikingly, they modify the details of cosmological nucleosynthesis, an effect 
we investigate in depth in this paper. PB Hs evapo rating at later times could also have important astrophysical effects, 
such as helping to reionize the Universe EH Hal. 

Even if PBHs had none of these effects, it is still important to study them because each one is associated with an 
interesting upper limit on the fraction of the mass of the Universe which can have gone into PBHs on some mass 
scale M . This fraction is epoch dependent but its value at the formation epoch of the PBHs is of great cosmological 
significance, this quantity being denoted by /3(M) . The current density parameter fipBH (in units of the critical 
density) associated with unevaporated PBHs which form at a redshift z or time t is roug hly related to /3 by 0] 

'tV 1/2 M ^ 1/2 



ftpBH^Wl + ^-lO 6 /? (^— J ~10 18 /3 f I p r -J (M>10 15 g), (1.4) 

where fi r ~ 10~ 4 is the density parameter of the cosmic microwave background (CMB) and we have used Eq. (jl.ip . We 
have also neglected the effect of entropy production after PBH formation. A more precise formula will be given later. 
The (1 + z) factor arises because the radiation density scales as (1 + z) 4 , whereas the PBH density scales as (1 + z) 3 . 
Any limit on SIpbh therefore places a constraint on /3(M) . For example, the 7-ray limit implies /3(10 15 g) < 10~ 26 
and this is one of the strongest constraints on j3 over all mass ranges. Another immediate constraint for PBHs with 
M > 10 15 g comes from requiring ^pbh to be less than the total matter density. As discussed in this paper, there 
are also many constraints on [3(M) for PBHs which have already evaporated, although the parameter f2pBH must 
then be reinterpreted since they no longer contribute to the cosmological density. Note that Eq. (|1.4j) assumes that 
the PBHs form in the radiation-dominated era, in which case /3 is necessarily small. Indeed, /3 could only reach 1 
for M ~ 10 17 Mq , which corresponds to the horizon mass at matter-radiation equality, and such black holes could 
hardly be regarded as primordial anyway. 

The constraints on (3(M) were first brought together by Novikov et al. 139]. An updated version of the constraints 
was later provided by Carr et al. [8J] , although this contained some errors which were corrected by Green and Liddlc 
[l30j |. Subsequently, the (3(M) diagram has frequently been updated and improved as the rele vant effects have been 
studied in greater detail. The most recent version of this diagram comes from Josan et al. [l40j |. The important 
qualitative point is that the value of (3(M) must be tiny over almost every mass range, even if the PBH density is 
large today, so any cosmological model which would entail an appreciable fraction of the Universe going into PBHs is 
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immediately excluded. F or exampl e, this places strong constraints on the amplitude of the density inhom ogeneitie s 
in the early Universe [3, Ill3l Il4l| and on the deviations of such inhomogeneities from Gaussianity 
One can also infer indirect limits on the spectral index of the primordial density fluctuations [H (U Il30l 1461 - 148j | 
and constrain the reheating process which follows inflation [35| |36|. In a less con ventional context, one can cons train 



cosmological models involving a time-varying gravitational "constant or extra dimensions 

pMIsiJ- We 



do not explore such scenarios in the present paper but it should be stressed that the form of the /3(M) constraints 
itself changes in this situation. 

One of the main purposes of the present paper is to update the limits associated with big bang nucleosynthesis 
(BBN) and the extragalactic 7-ray background (EGB) because these turn out to be the dominant ones over the mass 
range 10 9 -10 17 g. Both of these limits have been a subject of long-standing interest but it is necessary to reassess 
them in the light of recent observational and theoretical developments. On the observational front, there are new 
data on the light-element abundances, the neutron lifetime, and the 7-ray background. On the theoretical front, our 
view of the physics of strong interactions has changed. It used to be assumed that hadrons arc emitted from PBHs in 
the form of nucleons or mesons. However, according to the modern view of quantum chromodynamics (QCD), PBH s 
must produce quark and gluon jets which then fragment into hadrons over the QCD distance, Aqq D ~ 1CP 13 cm [l6l| . 
These hadrons then decay into the stable elementary particles of the standard model. 

The consequences of quark and gluon emission for the EGB limit were first studied by MacGibbon and Carr jl62| | 
but the analysis can now be refined because of convergence onto a standard (ACDM) cosmological model and because 
of improved calculations of the production of photons through jet decays. As regards the consequences of quark and 
gluon emission for BBN, most of the hadrons created decay almost instantaneously compared to the nucleosynthesis 
time scale, but long -lived ones remain in the ambient medium long enough to leave an observable signature. Kohri 
and Yokoyama [1631 ] first analyzed this effect but only in the context of the relatively low-mass PBHs evaporating in 
the early stages of BBN. In this paper we incorporate the effects of heavier PBHs, which evaporate after BBN, and 
infer constraints imposed by the dissociation and overproduction of synthesized light elements by both hadrons and 
high-energy photons. Photons from unstable mesons or hadrons and bremsstrahlung emission from quarks and gluons 
are also appropriately incorporated. We thus calculate for the first time the PBH constraints over the full relevant 



mass range. The code employed in our calculations is similar to the one used by Kawasaki et al. 164j in studying 
the effects of decaying particles on BBN. 

The plan of the paper is as follows. Section [H] describes the background equations and defines various quantities 
related to PBHs. Section [Til] reviews black hole evaporation and the effects of quark-gluon emission. Section IIVI 
then discusses the constraints deriving from cosmological nucleosynthesis effects, while Sec. [V] discusses the ones 
associated with the photon background. Section IVT1 combines both constraints in a single j3(M) diagram for the mass 
range 10 9 -10 17 g and then discusses some other (mainly less stringent) constraints in this mass range. Section IVIII 
summarizes the most important limits in mass ranges associated with larger nonevaporating PBHs. Section IVIIII 
collects all the constraints together into a single "master" (3(M) diagram and draws some general conclusions. It 
should be stressed that Sees. IVII and IVIII include quite a lot of review of previous work but it is useful to bring 
all the results together and we have elucidated earlier work where appropriate. Throughout most of this paper we 
assume that the PBHs have a monochromatic mass function, but allowing even a small range of masses around 10 15 g 
would have interesting observational consequences, especially for the EGB limits. However, this discussion is rather 
technical, so it is relegated to an Appendix. 



II. DEFINITIONS AND DESCRIPTION OF BACKGROUND COSMOLOGY 



In this section, we present some relevant definitions and background equations. We assume that the standard 
ACDM model applies, with the age of the Universe being tp = 13.7 Gyr, the Hubble parameter being h — 0.72 and 
the time of photon decoupling being tdec = 380 kyr |l65l Il66| . Throughout the paper we put c = H = kn = 1. The 
Friedmann equation in the radiation era is 

2 SttG 4n 3 G 4 
where g* counts the number of relativistic degrees of freedom. This can be integrated to give 



where and T are normalized to their values at the start of the BBN epoch. Since we are only considering PBHs 
which form during the radiation era (the ones generated before inflation being diluted to negligible density), the initial 
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PBH mass M is related to the "standard" particle horizon mass Mph (which is not the actual particle horizon mass 
in the inflationary case) by 



M = 7 M PH 



4tt 



jpH~ 3 « 2.03 x 10 5 7 



Is 



Mr, 



(2.3) 



Here 7 is a numerical factor which depends on the details of gravitational collapse. A simple analytical calculation 
suggests that it is around (l/\/3) 3 ~ 0.2 during the radiation era [H, although the first hydrodynamical calculations 
gave a somewhat smaller value [5[ . The favored value has subsequently fluct uated as people have performed more so- 



phisticated computations but now seems to have returned to the original one [167] . However , as mentioned earlier, the 
effect of critical phenomena could in principle reduce the value of 7 , possibly down to 10~ 4 [l68| . as could a reduction 
in the pressure |44l446l |. On the other hand, if the overdensity from which the PBH forms is "noncompensated" (i.e. 
not surrounded by a corresponding underdensity) , subsequent accretion could generate an eventual PBH mass well 
above the formation mass, leading to an "effective" value of 7 much larger than 1 Q. In view of the uncertainties, 
we will not specify the value of 7 in what follows. 

Throughout this paper we assume that the PBHs have a monochromatic mass function, in the sense that they all 
have the same mass M . This simplifies the analysis considerably and is justified providing we only require limits on 
the PBH abundance at particular values of M . Assuming adiabatic cosmic expansion after PBH formation, the ratio 
of the PBH number density to the entropy density, tt-pbh/s , is conserved. Using the relation p = 3sT/4, the fraction 
of the Universe's mass in PBHs at their formation time is then related to their number density ?ipbh(£) during the 
radiation era by 



p{M) 



_ PPBH^i) _ MnpBli(ti 



p(ti 



7.98 x 10" 29 7 



-1/2 



f 9*i \ 
V 106.75/ 



4M np BH (t) 
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s(t) 
M 



/ ^pbh(*0/ 
V 1 Gpc" 3 



(2.4) 



where the subscript "i" indicates values at the epoch of PBH formation and we have assumed s = 8.55 x 10 85 Gpc~ 3 
today. g*i is now normalized to the value of g* at around 10 s since it does not increase much before that in the 
standard model and that is the period in which most PBHs are likely to form. The current density parameter for 
PBHs which have not yet evaporated is given by 



n 



PBH — 
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(2.5) 



which is a more precise form of Eq. (jl.4l) . The dependence on the Hubble parameter is included here but we will 
generally assume h — 0.72. An immediate constraint on j3(M) comes from the limit on the CDM density parameter, 
CDM h 2 = 0.110 ± 0.006 with h = 0.72, so the 3ct upper limit is ^ PB h < ^cdm < 0.25 [UlIIll- This implies 



/3(M) < 2.04 x 10~ 18 7~ 1/2 



^CDM 
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(M > 10 ii5 g) 



(2.6) 



This constraint, which applies only for PBHs which have not evaporated yet, agrees with the one given by Blais et 
al. [20| when the parameters are appropriately renormalized. Note that the dependences on 7 and <?* in Eq. (|2 .4[) 
and subsequent equations arise through the relationship between M and 7} implied by Eqs. (|2.2j) and (|2.3[) . Since j3 



always appears in combination with 7 ' g 



1/2 -1/4 



it is convenient to define a new parameter 



P{M) = 7 1/2 



\ 106.75/ v ' 



(2.7) 



where g^ can be specified very precisely but 7 is rather uncertain. Most of the constraints discussed in this paper 
will therefore be expressed in terms of /3' rather than /3 . 

One can also represent the PBH constraints in terms of the ratio of the PBH density to radiation density at the 
evaporation epoch. This fraction is larger than the formation fraction by a redshift factor and given by 



a{M) = P{M) 



1 



1 -f - -^cvap 



/3(M) max 



M 



M V\ J \ Mpi Me. 



(2. 



where M cq is the mass of a PBH evaporating at the time of matter-radiation equality. The first and second expressions 
apply for PBHs which evaporate before and after then, respectively. A more precise expression has been derived by 
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Green and Liddlc [130]. They also distinguish between the fraction of the radiation density (/3) and the fraction of 
the total density (3/(1 + (3) going into the PBHs, but the difference is negligible so long as (3 is small. In any case, we 
just refer to f3(M) in this paper. 

Note that the relationship between f3 and fipBH must be modified if the Universe ever deviates from the standard 
radiation-dominated behavior. For example, if there is a dustlike stage for some extended early period t\ < t < ti , 
then one must include an extra factor (^Ai) 1 / 6 on the right- hand-side of Eq. ()2.6|1 |46| . This is also the period 
in which 7 is likely to be sma ll. The expression for (3(M) may also be modified in some m ass r anges if there is a 
second infla tiona ry phase [130l | or if there is a period when the gravitational constant varies [150] or there are extra 
dimensions jl60l ]. 

In principle, the analysis can be extended to the (more realistic) case in which the PBHs have an extended mass 
spectrum by interpreting f3(M) as the fractional density of PBHs at formation over a logarithmic mass interval. More 
precisely, we can write 

»pbh(AM) = dlnM , (2.9) 

where NpBn(M,t) is the cumulative number density of PBHs integrated over the mass range from to M . For a 
monochromatic mass function, in which all the PBHs have mass m, we then have Npbh(M, t) — tipbh(£) S(M — m) , 
where 9 is the Heaviside function. However, this more general analysis is only necessary if we wish to infer constraints 
which depend on the PBH mass spectrum itself. Although a precisely monochromatic mass spectrum is clearly not 
physically realistic, one would only expect the mass function to be very extended if the PBHs formed from exactly 
scale- invariant density fluctuations [3]; it is not clear whether this ever applies and it is certainly not expected in an 
inflationary scenario. However, it is important to stress that there are some circumstances in which the spectrum 
would be quite extended and this has the important implication that the constraint on one mass scale may also imply 
a constraint on neighboring scales. For example, we have mentioned that the monochromatic assumption may fail 
badly if PBHs form through critical collapse and the way in which this modifies the form of (3(M) has been discussed 
by Yokoyama 79] . Another important point is that if the PBHs with M ss M* have a spread of masses AM <~ M* , 
one would expect the evaporation process itself to lead to a residual spectrum with npBn(M, to) oc M 3 for M <C M* 
[Ifl ; the implications of this are discussed in the Appendix. 



III. EVAPORATION OF PRIMORDIAL BLACK HOLES 



A. Lifetime 

Here we briefly summarize some basic results of PBH evaporation. As first shown by Hawking [H], ISfl , a black hole 
with mass M = M\q x 10 10 g emits thermal radiation with temperature 

T BH = ^« 1 .06M-TeV. (3.1) 

This assumes that the hole has no charge or angular momentum, which is reasonable since ch arge and angular 
momentum will also be lost through quantum emission on a shorter time scale than the mass [l69Ml7l| . More 
precisely, such a black hole emits particles with energy between E and E + dE at a rate 

dl *s = =- ~~e7t V-TT27 ' ( 3 - 2 ) 

2tt e B / T »H - (-l) 2s 

where s is the spi n of the particle and T s is its dimensionless absorption coefficient, whose functional shape is found 
in Refs. [I69l4l71 |. It is related to the absorption cross section a s (M, E) by 



r s (M,E) = E2 ^ M > E K (3.3) 

In the high-energy limit, E ^> Tbh , o~s approaches the geometric optics limit er g = 27irG 2 M 2 . However, it is a 
more complicated function of E and M at low energies and depends on the spin. The average energies of the emitted 
neutrino, electron, and photon are E v = 4.22 Tbh, E e = 4 .18 T bh , and E 1 = 5.71 Tbh, respectively. The peak 
energies of the flux and power are within 7 % of these values |l6l| . 
The mass-loss rate of an evaporating black hole can be expressed as 

= -5.34 x 10~ 5 f{M) Afr 2 s" 1 . (3.4) 
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Here f(M) is a measure of the number of emitted particle species, normalized to unity for a black hole with M ^> 10 17 g, 
this emitting only particles which are (effectively) massless: photons, three generations of neutrinos and antineutrinos, 
and gravitons. The contribution of each relativistic degree of freedom to f(M) depends on the spin s [1721 ]: 

/ s=0 = 0.267, f s=1 = 0.060, f a=3/2 = 0.020, f a=2 = 0.007, 

f s =i/2 = 0.147 (neutral) , f s=1/2 = 0.142 (charge ±e) . '"" >l 

Holes in the mass range 10 15 g < M < 10 17 g emit electrons but not muons, while those in the range 10 14 g < M < 
10 15 g also emit muons, which subsequently decay into electrons and neutrinos. The latter range is relevant for the 
PBHs which are completing their evaporation at the present epoch. 

Once M falls to roughly 10 14 g, a black hole can also begin to emit hadrons. However, hadrons are composite 
particles made up of quarks held together by gluons. For temperatures exceeding the QCD confinement scale, Aqcd = 
250-300 MeV, one would expect these fundamental particles to be emitted rather than composit e par ticles. Only pions 
would be light enough to be emitted below Aqcd • Above this temperature, it has been argued that the particles 
radiated can be regarded as asymptotically free, leading to the emission of quarks and gluons. Since there are 12 
quark degrees of freedom per flavor and 16 gluon degrees of freedom, one would expect the emission rate (i.e., the 
value of /) to increase suddenly once the QCD temperature is reached. If one includes just u , d , and s quarks and 
gluons, Eq. (|3 . 5[) implies that their contribution to / is 3 x 12 x 0.14 + 16 x 0.06 « 6, compared to the pre-QCD 
value of about 2. Thus the value of / roughly quadruples, although there will be a further increase in / at somewhat 
higher temperatures due to the emission of the heavier quarks. After their emission, quarks and gluons fragment 
into further quarks and gluons until they cluster into the observable hadrons when they have travelled a distance 
Aq£ D ~ 10~ 13 cm. This is much larger than the size of the hole, so gravitational effects can be neglected. Note that 
the peak energy of the PBH emission reaches Aqcd once M falls to a critical mass M q s» 2 x 10 14 g, although the 
large increase in f(M) only occurs at a somewhat smaller mass. 

If we sum up the contributions from all the particles in the standard mod el up to 1 TeV, corresponding to M\q ~ 1, 
this gives f(M) — 15.35. Integrating the mass-loss rate over the time |l72| |. we find a lifetime 

T ~407 ( ffi-H 1 MfoS. (3.6) 



, 15.35, 

This can be inverted to give the mass of a PBH evaporating at time t after the big bang 



M « 1.35 xW^Y /3 m 1/3 g. (3.7) 



V 15.35 y Vis 

A PBH evaporating before the end of BBN at t ~ 10 3 s therefore has M < 10 10 g and Tbh 1 TeV. The critical mass 
for which r equals the age of the Universe is denoted by M* . If one assumes the currently favored age of 13.7 Gyr, 
one finds 

M, « 1.02 x 10 15 1 g « 5.1 x 10 14 g, (3.8) 

where /* is the value of / associated with the temperature Tbh(M*) • The indicated value of M* corresponds to 
/* = 1.9 and Tbh = 21 MeV. At this temperature muons and some pions are emitted, so the value of /* accounts for 
this. Although QCD effects are initially small for PBHs with M — M* , only contributing a few percent, it should 
be noted that they will become important once M falls to M q w 0.4 M* since the peak energy becomes comparable 
to Aqcd then. As discussed below, this means that an appreciable fraction of the time-integrated emission from the 
PBHs evaporating at the present epoch goes into quark and gluon jet products. 

It should be stressed that the above analysis is not exactly correct because the value of f(M) in Eq. (|3.7p should really 
be the weighte d ave rage of f(M) over the lifetime of the black hol e. T his accords with the more precise calculation 
of MacGibbon [l72l |. In fact, the update of the calculation in Ref. [173], using the modern value r = 13.7 Gyr, gives 
M* = 5.00 x 10 14 g, which is slightly smaller than our value. However, we will continue to use Eq. ()3.7|1 because we 
need the rough scaling with f(M) later in the paper. The weighted average of f(M) is well approximated by f(M) 
unless one happens to be close to a particle mass threshold, so the error involved in our approximation is usually 
small. For example, since the lifetime of a black hole of mass 0.4M, is roughly 0.25 x (0.4) 3 = 0.016 that of an 
black hole, one expects the value of M* to be overestimated by a few percent. This explains the small difference 
between our own and MacGibbon's more precise calculation. 
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B. Particle spectra 



Particles injected from PBHs have two components: the primary component, which is the direct Hawking emission, 
and the secondary component, which comes from the decay of the hadrons produced by fragmentation of primary 
quarks and gluons, and by the decay of gauge bosons. For example, the total photon spectrum emitted by a low-mass 
PBH can be written as 

diV^ dAT ri diV" cc 

^(E„M) = —j-^—(E 1 ,M) + -^L-{E^M) , (3.9) 

with similar expressions for other particles. In order to treat QCD fragmentation, we use the PYTHIA code, a Monte 
Carlo event generator jl74j | . This code has been constructed to fit hadron fragmentation in the high-energy regime, in 
particular reproducing the latest experimental results for center-of-mass energies below 200 GeV. We use this package 
for even larger center-of-mass energies by assuming that standard model physics can still be applied. It also gives 
the energy spectrum of created particles such as nucleons and pions, and photons from unstable mesons/hadrons and 
bremsstrahlung during the hadron fragmentation are now appropriately incorporated. However, the bremsstrahlung 
from secondary charged leptons and long-lived charged hadrons is not considered since it is subdominant. 



The spectrum of secondary photons generated by PBHs was first analyzed by MacGibbon and Webber [161] , It 
is peaked around E y ~ m T o /2 w 68 MeV, independent of the Hawking temperature, because it is dominated by the 
27-decay of soft neutral pions which are practically at rest. As E 1 increases, it then falls as E~ x before cutting off as 
E^ exp(— -E 7 /Tbh) above Tbh • The peak flux can be expressed as 

— ^(£ 7 =m ff o/2)~2— ^(E„ =m„o)~2 V B^o(E 7 m„«) ^-( Ei = E), (3.10) 

where B q g ^ n o (£jot, E^o) is the fraction of the jet energy Ej ct going into neutral pions of energy E„o . This is of 
order 0.1 and fairly independent of jet energy. The second approximation in Eq. (|3 . 10[) uses the fact that most of the 
primary particles have of order of the average energy E « 4.4 Tbh • From Eq. (|3.2|) . we then find that 

(with units of /zT 1 ) for any value of s in the geometric optics limit. Thus the energy dependence of Eq. (|3.10|) comes 
entirely from the factor E ss 4.4 Tbh and is proportional to the Hawking temperature. A similar analysis applies for 
the secondary production of protons and antiprotons, this also scaling as the black hole temperature. 

The emission rates of primary and secondary photons for four typical temperatures are shown in Fig. [TJ Thes e 
results might be compared with those obtained by HERWIG |l75| . the code used by MacGibbon and Webber [l6lj . 
Although HERWIG and PYTHIA have been continually updated to match almost 20 years of high-energy accelerator 
runs, it should be stressed that there have been few new e + e~ accelerator experiments near Aqcd ■ Most new e + e~ 
experiments have been at the highest accelerator energies, while recent experiments examining behavior near Aqcd 
have involved heavy ion (i.e. quark-gluon plasma) rather than e + e~ collisions, the latter being more closely analogous 
to the black hole situation. Hence there has been little recent advance in our understanding of processes near Aqcd 
and any accelerator code differences probably come from different choices for u and d quark masses and the value of 
Aqcd • These parameters can be changed arbitrarily in those codes anyway. 

If one integrates the instantaneous spectrum dN 1 /AE 1 over time for M < , one obtains the more observationally 
significant quantity dA 7 /d£' 7 . At most energies this is dominated by the emission associated with the initial PBH 
mass M , since the mass only changes in the relatively short period towards the end of evaporation. However, PBHs 
with M < M* also have a high-energy tail for E 1 3> Tbh (M) due to their final phase of evaporation. (For observations 
at the present epoch, this energy must be reduced by a factor l + z cvap (M) where Zevan(M) is the evaporation redshift.) 
The precise form of this high-energy tail has been investigated by MacGibbon [l72| | and it is discussed qualitatively 
in the Appendix, so we just summarize the conclusion here. For M > M q , one obtains 

d/^ [E* (£ 7 <m-i), 

dE^lE-z (Af 1 < E. t < Aqcd) , 1 ' 

where the lower energy part comes from the initial Ray leigh- Jeans contribution. The form of the tail changes for 
.Ely > Aqcd due to secondary emission, involving an intermediate E~ x regime, but it still falls off as E~ 3 at very high 
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FIG. 1. Instantaneous emission rate of photons for four typical black hole temperatures. For each temperature, the curve with 
the peak to the right (left) represents the primary (secondary) component and the thick curve denotes their sum. 



energies. For M < M q , the spectrum is always dominated by the secondary photons. As M increases towards M* , 

the primary photons generate a low-energy redshifted tail which scales as El/ 2 |l72| . 

It should be noted that ratio of the time-integrated secondary flux to primary flux drops rapidly once M goes above 
M* . This can be understood as follows. We have already seen that a black hole with M = M» will emit quarks 
efficiently once its mass gets down to M q ss 0.4 and this corresponds to an appreciable fraction of its original mass. 
On the other hand, a simple calculation shows that a PBH with somewhat larger initial mass, M = (1 + fj,) M* , and 
hence a slightly longer lifetime, r w (1 + /i) 3 to , will today have a mass 



M(t ) « (3 M ) 1/3 (1 + M + M73) i/3 M, 



1/3 



(3.13) 



providing /i is not too small (as explained below). Here we have assumed f(M) ~ /* , which should be a good 
approximation for m > M q since the value of / only changes slowly above the QCD threshold. For /^>1, Eq. (|3.13|) 
just gives m w fi w M , as expected. For the second term can be neglected and m exceeds M q providing 

(Sfj,) 1 / 3 > 0.4, which corresponds to \i > 0.02. Equation (|3. 13[) must be modified for fj, < 0.02 since the value of / 
increases by a factor of around 4 as M falls below M q . If we assume that it jumps discontinuously from /» to a /„ at 
M q , then the mass falls to M q at a time 



1- (0.4) 3 + 3 fi 
1- (0.4) 3 (1 -a- 1 ) 



to 



(3.14) 



and this leads to a current mass 



m^(3«/i) 1/3 M» (/i<0.02). 



(3.15) 



This differs from the expression given by Eq. (|3.13[) because of the a factor but it goes to zero at fx = 0, as expected. 
The denominator in the expression for t q arises because the lifetime of a black hole of mass close to M* is slightly 
increased by the change of / . 

The fact that m falls below M q only when fx is less than the tiny value 0.02 means that the fraction of the black 
hole mass going into secondaries falls off sharply above M* . Indeed, this fraction must be around 0.4 (1 — m/M q ) 
for fi < 0.02, since the unevaporated mass m does not contribute, and it must be exponentially reduced by the 
Wien blackbody factor oc exp(— \ m/Mq) for fi > 0.02, where x is a parameter which depends only weakly on ji. A 
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transition occurs between the /i < 0.02 and fi > 0.02 regimes. The ratio of the secondary to primary peak energies 
and the ratio of the time-integrated fluxes are shown in Fig. [5J they are qualitatively well explained by the above 
arguments. Note that the forms of the curves at around log 10 (M/g) = 13.4 are associated with interactions of photons 
with the background Universe at high redshifts. These are discussed later and affect the secondary particles more 
than the primary ones. 
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FIG. 2. Ratios of the time-integrated secondary to primary peak energies (solid red) and fluxes (dashed green). 



C. Photosphere effects 



There has been some dispute in the literature about the interactions between emitted particles around an evaporating 
black ho le. T he usual assumption |l76| that there is no interaction between emitted particles has been disputed by 
Heckler [177| . who claims that QED interactions could produce an optically thick photosphere once the black hole 
temperature exceeds Tbh = 45 GeV. He ha s proposed that a similar effect may operate at an even lower temperature, 
Tbh ~ 200 MeV, due to QCD effects Il78ll. Variants of these mo dels and t heir astrophysical implications have been 
studied by Cline et al. |l79j |. Kapusta jl80l | and Daghigh et al. [I8ll4l83| . However, MacGibbon et al. |l73| have 
analyzed these claims and identified a number of physical and geometrical effects which invalidate them. There are 
two key problems. First, the particles must be causally connected in order to interact and only a negligible fraction 
of the emitted charged particles satisfies this constraint. The standard cross sections are further reduced because the 
particles are created at a finite time and do not go back to the infinite past. Second, a scattered particle requires 
a minimum distance to complete each bremsstrahlung interaction, with the consequen ce th at there is unlikely to be 
more than one complete bremsstrahlung interac tion per particle near the black hole |l84j . This relates to what is 
termed the Landau-Pomeranchuk-Migdal effect 185l4l87j | . 



MacGibbon et al. conclude that the emitted particles do not interact sufficiently to form a QED photosphere 
and, using analogous QCD arguments, that the conditions for QCD photosphere formation could only be temporarily 
satisfied (if at all) when the black hole temperature is of order the QCD confinement scale, Aqcd = 250-300 MeV. 
Even in this case, the strong damping of the Hawking production of QCD particles and the limited available energy 
per particle around this threshold should suffice to suppress it. In any case, no QCD photosphere per sists once 
the black hole temperature climbs above Aqcd • They also consider the suggestions of Belyanin et al. 188] that 
plasma interactions between emitted particles could form a photosphere but conclude that this too is implausible. A 
final possibility is that some form of string photosphere might arise due to "stretched horizon" effects, as discussed 
in Ref. [l89j |. but this would only apply at energies too high to be of interest here. Throughout this paper, we 
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therefore usually adopt the standard assumption that no photosphere ever forms, although we do return to this issue 
in discussing PBH explosions. 

It must be stressed that there may well be interactions between emitted particles in the background Universe even 
if there is no interaction near the black hole, so there is a sense in which one may have a cosmological photosphere. 
Although there is no cosmological QCD photosphere in the period after 10~ 5 s relevant to BBN, there will be a 
cosmological QED photosphere before around 10 s since background electrons and positrons will then be as numerous 
as photons. However, the situation is different in the cosmological photosphere context because the emitted particles 
are not interacting with each other. Also, in the black hole photosphere context, the particle density is much higher 
near the hole and falls off with the square of distance until it reaches the background cosmological value. 

IV. CONSTRAINTS ON PBHS FROM BIG BANG NUCLEOSYNTHESIS 

A. Historical overview 

In this section we will be concerned with the effects of PBH evaporations on BBN, presenting for the first time the 
full analysis of the associated constraints on the PBH mass spectrum. There is also an indirect limit on much larger 
mass scales associated with PBH formation because the density inhomogeneities at around neutron-proton freeze-out 
would modify BBN if they were large enough to generate PBHs. This corresponds to mass scales around 10 5 Mq but 
the associated limit is not definitive because it assumes a particular formation mechanism for the PBHs. 

Since the effects of PBH evaporation on BBN have been a subject of long-standing interest, we start with a brief 
review of previous work and then explain why it is necessary to reassess the c onstr aints. All the limits will be expressed 
in terms of the parameter (3'(M) defined by Eq. (I2.7[) . Vainer and Nasel'skii |190l ] studied the effects of the injection of 
high-energy neutrinos and antineutrinos from PBHs. This changes the epoch at which the weak interactions freeze out 
and thereby the neutron-to-proton ratio at the onset of nucleosynthesis. This results in an increase of He production, 
so by demanding that the primordial abundance satisfies Y p < 0.33, the most conservative constraint available at the 
time, they concluded that the ratio of the PBH density to the baryon density at BBN should be smaller than 10-10 4 
for M — 10 9 -3 x 10 11 g. Equation (|2.5I) implies that this corresponds to a limit 

P'(M) < 3 x (10~ 18 -10~ 15 ) M\l 2 ( M = lo9 - 3 x 1Ql1 g) • (4- 1 ) 

These authors also studied the effects of entropy generation by PBHs but obtained only the modest constraint that 
their density should be less than that of the backg roun d radiation at evaporation. More detailed numerical analysis 
of this effect was carried out by Miyama and Sato |l9lj , who calculated the entropy production from the PBHs with 
AI = 10 9 -10 13 g which evaporated during or after BBN. If these PBHs were too abundant, the baryon-to-entropy 
ratio at nucleosynthesis would be increased, which would result in overproduction of 4 He and underproduction of D. 
They demanded that the primordial mass fractions of these elements satisfy Y p < 0.29 and D p > 1 x 10~ 5 . One can 
approximate their constraint with a single power law: 

P'(M ) < 10~ 15 M w 5/2 (AI = 10 9 -10 13 g) . (4.2) 

Zel'dovich et al. [l92j studied the effect of PBH emission of high-energy nucleons and antinucleons. They claimed 
that this would increase the deuterium abundance due to capture of free neutrons by protons and spallation of 4 He. 
They inferred the upper limits: 



P'(AI) < 



6 x 10" 


-18 


M W 112 


(M = 


10 9 -10 10 g) , 


6 x 10" 


-22 




(M = 


10 10 -5 x 10 10 g) , 


3 x 10- 


-23 


mT 


(AI = 


5 x 10 10 -5 x 10 n g), 


3 x 10- 


-21 


M w 1/2 


(AI = 


5 x 10 n -10 13 g) . 



(4.3) 



Vainer et al. [193l | then performed numerical studies of neutron-antineutron injection using a nucleosynthesis network. 
They found that the spallation of 4 He resulted in extra D production. They took the neutron lifetime to be 918 (±14) s 
and used the observational results Y p — 0.29 ±0.04 and D p = 5 x 10~ 5 . On the assumption that PBHs were produced 
from scale-invariant density fluctuations, rather than ha ying a monochromatic spectrum, they obtained the constraint 
f3' < 10 -26 for the relevant mass range. Later Lindley [l94j | considered the photodissociation of deuterons produced 
in nucleosynthesis by photons from evaporating PBHs with M > 10 10 g. He found the constraint 

P'(M) < 3 x 10~ 20 Mi /2 (M > 10 10 g) (4.4) 
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and con cluded that photodestruction was comparable to the extra production of deuterons discussed by Zel'dovich et 
al. [l92j . All the above limits are shown in Fig. 1 of Ref. Q . 

Observational data on both the light-element abundances and the neutron lifetime have changed considerably since 
these classic papers were written. Much more significant, however, are developments in our understanding of the 
fragmentation of quark- ant iquark pairs or gluons emitted from a PBH into hadrons within the QCD distance. Most 
of the hadrons created decay almost instantaneously compared to the time scale of nucleosynthesis, but long-lived ones 
(such as pions, kaons, and nucleons) remain long enough in the amb ient medium to leave an observable signature on 
BBN. These effects were first discussed by Kohri and Yokoyama |163l ] but only in the context of the relatively low-mass 
PBHs evaporating in the early stages of BBN. In the rest of this section, we extend and update their calculations. 



B. Particle injection 



As mentioned above, in the Hawking temperature range of interest here, the relevant high-energy hadrons are 
not primary particles but decay products from hadron fragmentation. At first these hadrons scatter off background 
photons and electrons because these particles are much more abundant than background nucleons. We incorporate 
thermalization through the electromagnetic interactions described in Appendix B of Ref. [l64f . Specifically, for hadron 
species Hi , we take into account Coulomb scattering (Hi + — > Hi + e ), Compton scattering (Hi + 74^ + 7), 
the Bethe-Heitler process (iZj+7 — > Hi + e + + e~) and the photo-pion process (i/i+7 — > H^ + tt). Coulomb scattering 
dominates in the temperature regime of interest and thermalizes relativistic charged particles on a timescale 

r ch ~ 3 x KT 11 ( E \ ( - ^ s, (4.5) 

\1GcVJ V0.02McVy ' V ; 

where E is the kinetic energy of the charged meson [l95j |. This equals the hadron interaction timescale at a cos- 
mological temperature T ~ 0.04 m e sa 0.02 MeV, both times being much less than the cosmic expansion time then, 
so thermalization occurs above this temperature. For neutrons, Coulomb scattering occurs through their magnetic 
moment and the cross section is suppressed by a factor of T 2 jw? n , so they are thermalized only for T > 0.09 MeV. 

In the early stages of BBN, therefore, we expect emitted particles to be quickly thermalized and to have the 
usual kinetic equilibrium distributions before they interact with ambient nucleons. In addition, as we show later, it 
is reasonable to treat the emitted hadrons as being homogeneously distributed. In order to estimate the minimal 
time scale for strong interactions, we therefore use the thermally averaged cross sections (crv)^- i _^ N , for the strong 
interactions between hadrons and ambient nucleons N (i.e. protons or neutrons). For the hadron process N + Hi — > 
N' + ■ ■ ■ , the strong interaction rate is 

where 771 is the initial baryon-to-photon ratio (rib/n 7 at T > 10 MeV, where = n p + n n ), is the number density 
of the initial nucleon species, corresponding to a fraction /jv = tin /nh 1 ari d T„ is the neutrino temperature. Thus 
we only need to consider particles with lifetime longer than O(10~ 8 ) s. The corresponding mesons are 7r + , 7r~ , K + , 
K~ , and Kl , while the baryons are p , p , n , and n . In this work, we do not include effects induced by kaons because 
they are smaller than those induced by pions. 



C. Particle reactions 



High-energy particles emitted by PBHs would modify the standard BBN scenario in three different ways: (1) high- 
energy mesons and antinucleons induce extra interconversion between background protons and neutrons even after 
the weak interaction has frozen out in the background Universe; (2) high-energy hadrons dissociate light elements 
synthesized in BBN, thereby reducing 4 He and increasing D, T, 3 He, 6 Li, and 7 Li; (3) high-energy photons generated 
in th e cas cade further dissociate 4 He to inc reas e the abundance of lighter elements even more. We analyze (1) following 
Ref. [l63l | and (2) and (3) following Ref. 164 1 , which studied the effects of the hadronic decay of long-lived massive 



particles on BBN in the context of decaying dark matter scenarios. Since the presentation in those references are very 
detailed, we do not repeat them in this paper but only focus on the points specific to PBHs. 
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The Boltzmann equations corresponding to processes (l)-(3) are schematically written as 
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dt 
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dt 
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(N = p,n), 



(Ai = D, T, 3 He, 4 He, 6 Li, 7 Li) , 



(4.7) 



where, as argued above, the relevant species emitted by PBH evaporation are 7r + , n~ , p , p , n , and n , as well 
as photons and neutrinos. In these equations the subscript "SBBN" denotes the standard big bang nucleosynthesis 
contribution, which is followed using the KAWANO code 4.1 [l96| with updated nuclear cross sections; the subscript 
"conv" means interconversion between background neutrons and protons due to strong interactions; "hadron" repre- 
sents both production and destruction through hadronic showers; and "7" denotes production or destruction through 
photodissociation. 



1. Interconversion of protons and neutrons due to mesons 



Long-lived charged mesons are thermalized and scatter off the ambient nucleons with a thermally averaged cross 
section before they decay. The cross section of the exothermic process is inversely proportional to the velocity v 
when it is small and so we can use the threshold value. For an emitted (anti)nucleon, the efficiency of thermalization 
is greater. Those (anti)nucleons are efficiently stopped by electron scattering until T ~ 0.09 MeV for neutrons and 
T ~ 0.03 MeV for protons, after which we may treat a nucleon-antinucleon pair as a kind of ther mal m eson, which 
induces interconversions via nn + p — > n ■ ■ ■ and pp + n — > p + • • • . For details, see Sec. Ill of Ref. |l63j . 



2. Hadrodissociation 



For T < 0.09 MeV for neutrons and antineutrons or T < 0.02 MeV for protons and antiprotons, high-energy nucleons 
may scatter off the background nuclei produced by BBN to induce hadrodissociation before losing kinetic energy 
through electromagnetic processes. More precisely, we must specify the frequency of the hadronic reaction Hi + Aj —¥ 
Ak which occurs while the high-energy hadron is reducing its energy through the electromagnetic interactions. Here 
we consider only the background protons, ph g , and background helium nuclei, ab g , as target nuclei Aj because they 
are the dominant nuclear content at the relevant epoch. 

There are three main effects in this regime: (a) interconversion of background protons and neutrons through hadronic 
collisions (i.e. by energetic inelastic pp and pn scattering rather than by pp and pn discussed above); (b) dissociation 
of background ab g , created in BBN, into energetic debris such as n , p , D, T, and 3 He; and (c) production of heavier 
nuclei, such as 6 Li, 7 Li, and 7 Be, through collisions between these debris and background a^ g ■ For neutrons, we 
take into account the fact that they may decay before scattering off the background nuclei. We consider 10 reactions 
relevant to (a) and (b) and these are summarized in Table HI As regards (c), we incorporate the nonthermal production 
of 6 Li by energetic T or 3 He through T + — > 6 Li + n and 3 He + — >• 6 Li + p . Collisions of a^ g with high-energy 
4 He may result in the production of 6 Li, 7 Li, or 7 Be, although they do not affect the final constraint on the mass 
spectrum of PBHs. 

The necessary procedure is then to evaluate the number of each stable nuclei species produced per single PBH 
evaporation and to compare this with observations. Since the reaction rate depends on energy, it is important to 
calculate the energy spectrum of injected hadron species, as well as the number density. In our analysis the energy 
spect ra o f primary protons and neutrons are calculated by PYTHIA. For more details, see Sees. VII and VIII of 
Ref. PH. 



3. Photodissociation 



Even if high-energy quarks and gluons mainly produce hadronic particles initially, their kinetic energy is eventually 
converted into high -energy radiation through scattering. Hence it is important to incorporate photodissociation. 
Following Ref. [l97| . we consider the following processes in addition to the injection of particles from PBHs: (i) 
double-photon pair creation, 7 + 7b g — > e + + ; (ii) photon-photon scattering, 7 + 7b g — > 7 + 7 ; (iii) Compton 
scattering off background thermal electrons, 7 + — > 7 + e~ ; (iv) inverse Compton scattering off background 
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TABLE I. Hadronic reactions with p bg and a bg included in our calculation. 



Process 


i 


= n 


i = p 


(i,Pb g ,l) 


n + Pbg 


— > n + p 


P + Pbg -> P + P 


(*,Pb g ,2) 


n+Pbg- 


n 

> n + p + 7T 




P + P bg ->P + P + 7T 


(i,Pb g ,3) 


n + p bg ^ 


n + n + 7r T 


_l_ 

P + Pb g — ► P + « + 7T 


(j,Q bg ,4) 


n + a bg 


— 5- n + a 


p + a bg — > p + a 


(i (Tu M 


71 -t- C V i 
' T ! - t bg 


— > D + T 


n J- fV K — ^ D 4- 3 He 


(j,a bg ,6) 


n + a bg - 


-> p + n + T 


p + a bg — 5- 2p + T 


(j,a bg ,7) 


+ a bg 


n + 2D 


p + a bg -> P + 2D 


(j,Q bg ,8) 


n + ab s — 


► p + 2n + D 


p + Qb g — > 2p + n + D 


(j,a bg ,9) 


n + a bg 


-> 2p + 3n 


P + Qb g — >■ 3p + 2n 


(i,a b g,10) 


n + a bg — 


> n + a + 7T° 


p + abg — ► p + a + 7T° 



photons, e + 7b g — > e + 7 ; (v) pair-creation off background protons, 7 + pbg — > e + + e + p . The relevant terms 
in the Boltzmann equations are 

= y)/ d£ 7 (TA,^(£ 7 )/ 7 (£ 7 ) + VnA ] / d£ 7 a Aj ^4 ( (-E 7 )/ 7 (£ 7 ), (4.8) 

where ct^.^^ . is the cross section for the reaction A4 + j Aj + ■ ■ ■ with threshold energy £ 7 th) and / 7 (£? 7 ) is 
the photon energy distribution function, which is calculat ed ta king above processes (i) through (v) into account. 
The resulting spectrum has a cutoff at £™ ax ~ 771^/(22 T) [l98j |. above which photons lose energy through the pair 
creation (i). If £™ ax exceeds the threshold energies of light elements, high-energy photons can dissociate them and 
change their abundances. Since the binding energies of D and 4 He are Bd = 2.2 MeV and B a = 20MeV, they are 
photodissociated only after 10 4 s and 10 6 s, respectively. Specifically, we incorporate the following photodissociation 
reactions: D + 7 -> n+p; T + 7 n + D, p + 2n; 3 He + 7 -» p + D, 2p + n; 4 He + 7 -t p + T, n + 3 He, p + n + D. We 
also include photodissociation of 6 Li, 7 Li, and 7 Be. For more details, see Sec. IV of Ref. jl64| or Refs. [I97lll99l . l200j |. 



dn Ai 
dt 



D. Current status of observations of light elements 

The observed light-element abundances are listed in Table ILT1 and we now summarize the status of the observations. 



TABLE II. Observational constraints on abundances of light elements at 2a to be used in the present paper. 



Element Constraints Refs. 

4 He Yp = 0.2516 ± 0.0080 [201-2061 

D D/H < 5.16 x 10" 5 [207 . 2081 

3 He 3 He/D < 1.37 [209] 

6 Li 6 Li/ 7 Li < 0.302 [210 . 2111 



Helium. Synthesis of 4 He is the most important consequence of BBN and its abundance should in principle be 
used as the most fundamental test of nonstandard BBN. The primordial abundance of 4 He is inferred by using 
recombination lines from low-metallicity extragalactic HII regions which yield Hell/HII ratios. The primordial mass 
fraction, Y p , is obtained by extrapolating to th e O/ H — > or N/H — > limits. The recent tre nd is that Y p has a larger 
value than thought before. Olive and Skillman [203j | reanalyzed the seven HII regions of Ref. |20l| . using Hel emission 
lines to determine the t empe rature and including underlying stellar absorption, to obtain Y p = 0.2491 ± 0.00 91. 



202 to 



Fukugita and Kawasaki [2041 ] . on the other hand, used the OIII emission line in 33 HII regions in Ref. 
determine the temperature. Incorporating s tella r absorption, they found Y p = 0.250 ± 0.004 but Y p = 0.234 ± 0.004 
without stellar absorption. Peimbert et al. |205 | used new computational results of Hel emission lines to deter mine 
the temperature. From five HII regions of Ref . 201 |. they found Y p — 0.249 ±0.009. Most recently, Izotov et al. |206l ] 
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found Y p = 0.2516 ±0.0011 from 93 HII regions in HeBCD and 271 HII regions in SDSS DR5. In this paper, we adopt 
their mean value and the 2a error of Fukugita and Kawasaki [2041 ] : 

Y p = 0.2516 ±0.0080. (4.9) 

Deuterium and helium-3. The primordial deuterium-to-hydrogen abundance ratio D/H is measured in high redshift 
QSO absorption systems. Because such a Lyman-limit system is not expected to be contaminated by any Galactic 
or stellar chemical evolution, we can regard it as having primordial abundances. There have been several recent 
measurements with improved modelling of the continuum level, the Lyman-a forest, and the ve locity str ucture of the 
absorption systems. The weighted mean of these measurements, D/H = (2.82 ± 0.26) x 10 -5 [2071 . |208| . agrees well 
with the predi ction of standard BBN with the baryon-to-photon ratio inferred by observations of the CMB anisotropy 

[mdiiHil. Recentl y 

a new observatio n to wards Q0913+07, giving the same weighted value but a smaller error 
D/H = (2.82 ± 0.20) x 10" 5 was reported [H|. Among the data, the highest value is D/H = (3.98±g;|?) x 10~ 5 . In 
this paper we adopt this highest value to get a conservative 2a upper bound 

^ < 5.16 x 10~ 5 . (4.10) 
H 

As for 3 He, we use the upper limit on 3 He/D because thi s rat io is an increasing function of cosmic time. A recent 
measurement in protosolar clouds gives a 2a upper bound [209j ] 

3 He 

— < 1.37. (4.11) 

Lithium. The primordial value of 7 Li/H is measured in old Popula tion II ha lo stars. However, there is a discrepancy 
between the observational values and theoretical predictions in BBN |214Ll215 |. Therefore we do not adopt a constraint 
from 7 Li/H in this work. 6 Li is much rarer than 7 Li and its abundance is believed to increase with cosmic time 
through cosmic-ray spallation, so one can obtain an upper bound on its primordial abundance by observing low- 
metallicity regions. Asplund et al. [210l ] detected 6 Li in 9 out of 24 metal-poor halo dwarfs at more than 2a 
significance. In particular, they detected 6 Li in one very metal-poor star LP815-43 with [Fc/H] = —2.74 and inferred 
6 Li/ 7 Li — 0.046±0.022. Based on this result and a possible error (+0.106) discussed in Ref. [21 lj ] for stellar destruction 
processes, we take a conservative 2a upper bound 

6 Li 

TTj-r < 0.302. (4.12) 



E. Constraints on /3'(M) imposed by BBN 

We now report our results on the constraints on the PBH mass spectrum imposed by BBN. We have three param- 
eters: the initial baryon-to-photon ratio i]i , the PBH initial mass M or (equivalently) its lifetime r , and the initial 
PBH number density normalized to the entropy density, Ypbh = ^pbh/s . We start the BBN calculation at a cosmic 
temperature T = 100 MeV. The initial baryon-to-photon ratio needs to be set to an appropriate value, so we use the 
present one r\ — (6.225 ± 0.170) x 10~ 10 [166], after allowing for possible entropy production from PBH evaporations 
and photon heating due to e + e~ annihilations. From Eq. (|2.4p the value of Ypbh is related to the (normalized) initial 
mass fraction /3' by 

/3' = 5.4x 10 21 (^) 1/2 Ypbh- (4.13) 

In general j3' , r , and Ypbh in this relation would all depend on the PBH mass M , although we presuppose a 
monochromatic mass function in what follows. 

Figure |31 summarizes the results of our calculations. First, we note that PBHs with lifetime smaller than 10 _2 s 
are free from BBN constraints because they evaporate well before weak freeze-out and leave no trace. PBHs with 
M = 10 9 ~10 10 g and lifetime r = 10 _2 -10 2 s are constrained by the extra interconversion between protons and neutrons 
due to emitted mesons and antinucleons, which increases the n/p freeze-out ratio as well as the final He abundance. 
For t = 10 2 -10 7 s, corresponding to M = 10 10 -10 12 g, hadrodissociation processes become important and the debris 
deuterons and nonthermally produced 6 Li put strong constraints on /3(M) . Finally, for t = 10 7 -10 12 s, corresponding 
to M = 10 12 -10 13 g, energetic neutrons decay before inducing hadrodissociation. Instead, photodissociation processes 
are operative and the most stringent constraint comes from overproduction of 3 He or D. However, even these effects 
become insignificant after 10 12 s. 
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FIG. 3. Upper bounds on /3'(M) from BBN, with the broken lines giving the earlier limits in Ref. [16 



For comparison, we have also shown the much we aker constraint imposed by the entropy production from evapo- 
rating PBHs, as first analyzed by Miyama and Sato [l9ll |. The factor A in Fig. [3] is the ratio of the entropy density 
after and before PBH evaporations, which is calculate d num erically. We also show as dotted lines in Fig. [3] the limits 
obtained about ten years ago by Kohri and Yokoyama [l63j . It is noted that their helium limit on /3(M) was stronger 
than ours, while their deuterium limit was weaker. This is because the assumed helium abundance is now smaller, 
whereas hadrodissociation of helium produces more deuterium. 

It is important to note that the limit on f3(M) does not have a sudden cutoff above the mass Mi = 10 13 g associated 
with the PBHs which evaporate at t\ = 6.2 x 10 11 s. This is because even PBHs with mass M 2 larger than M\ will 
generate some emission before t\ . However, the fraction of their energy involved is only 



1 - 



1-1^ 

T2 



1/3 



Mi 
M2 



(4.14) 



^7/2 



for Mi < Mi , so Eq. (j2~4")l implies that the limit on /3(M 2 ) is weaker than that on /?(Mi) by a factor (M 2 /M] 
This explains the form of the upper cutoff above 10 13 g in Fig. [3l Indeed, this same argument implies that one never 
has a sudden upper cutoff in (3(M) for any limit which involves the tota l ene rgy emitted by the PBH. 
Finally we note that the neutron diffusion length at t = r is given by [216| 



d n ( T ) = 5.8 x 10 2 



T 



lMeV 



-5/2 



-1/2 



(4.15) 



where at is a tran sport cross section. From Eq. (I2.4p . the number of PBHs within the neutron diffusion volume is 
therefore given by [163] 



N d (r) = 7.7 x 10 46 M 3 d n (rf H{rf 



1.2 xlO 21 ^' (t = 1s), 

2.3 xlO 26 ^' (r=10 3 s) 



(4.16) 
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this scaling as Af~ 6 T~ 15 / 2 oc Af -6 t 15 / 4 oc t 7 / 4 . Comparison with Fig. [3] shows that this number always exceeds 1 
throughout BBN, so there are sufficiently many PBHs within the neutron diffusion length for the relevant parameter 
range for them to be observationally constrained. Therefore we are justified in assuming that the emitted hadrons 
form a homogeneous background. 



V. CONSTRAINTS ON PBHS FROM EXTRAGALACTIC PHOTON BACKGROUND 



A. Historical overview 



One of the earliest works that applied the theory of black hole evaporation to astrophysics was carried out by 
Page and Hawking [861 ] . who used the diffuse extragalactic 7-ray background observations to constrain the mean 
cosmological number density of PBHs which are completing their evaporation at the present epoch to be less than 
10 4 pc~ 3 , although their local density could be much larger than this if they are clustered inside the Galactic halo. 
This also correspon ds to an upper limit on their density parameter fipBH of around 10~ 8 . In fact, this limit was first 
noted by Chapline |217l ] and it was also studied by Carr [88| . though in less detail. The limit h as be en subsequently 
refined but without entailing large quantitative changes. For example, MacGibbon and Carr [l62j ] considered how 
it is modified by including quark and gluon emission effects and inferred fipBH < (7-6 ± 2.6) x 1CP 9 h~ 2 . (They 
gave the exponent of h as —1.95 ± 0.15 rather than —2, because of the M depe nden ce of /, but the difference is 
negligible.) Later they used EGRET observations between 30MeV and 120 GeV [218j ] to derive a slightly stronger 
limit fipBH < (5.1 ± 1.3) x 10~ 9 h~ 2 [Hi]. Using the modern value of h gives ft PB H < (9.8 ± 2.5) x 10~ 9 and this 
corresponds to the (h independent) constraint /3'(Af*) < 6 x 10~ 26 from Eq. (|2.5|) . They also inferred from the form 
of the 7-ray spectrum that PBHs do not provide the dominant contribution to the background. As discussed below, 
more recent EGRET and Fermi Large Area Telescope (LAT) data and more precise calculations of the PBH emission 
allow the limit on ^'(A/*) to be improved but the second conclusion remains true. This raises the point that one 
should in principle remove any known astrophysical sources of an extragalactic 7-ray background before c alcu l atin g 
the PBH constraints. For example, the contrib ution of blazars has recently been discussed by Inoue et al. 220, 221]. 
This is the strategy adopted by Barra u et al. [222j |. who thereby obtain a limit fipBH < 3.3 x 10~ 9 which is three 
times stronger than that in Ref. [219j . We do not attempt to make such a subtraction here but merely emphasize 
that our constraints on /?'(M) are therefore more conservative than necessary. 



B. Photon spectra 



The photon emission has a primary and secondary component, as indicated by Eq. (|3 . 9[) . The two spectra are 
calculated according to the prescription of Sec. IIIIBI One of our purposes here will be to determine the relative 
magnitude of these two components. This is sensitive to the PBH mass and affects the associated /3(M) limit. In 
order to determine the present background spectrum of particles generated by PBH evaporations, we must integrate 
over both the lifetime of the black holes and their mass spectrum, allowing for the fact that particles generated in 
earlier cosmological epochs will be redshifted in energy by now. This is quite a complicated calculation because the 
PBH mass is itself a function of time due to evaporation. However, for the purpose of obtaining limits on (3'(M) , we 
are assuming a monochromatic mass spectrum, which considerably simplifies the analysis. 

If the PBHs all have the same initial mass M , then the present-day photon flux is a superposition of the instanta- 
neous emissions from all previous epochs. If we approximate the number of emitted photons in the logarithmic energy 
bin AE 7 ~ E 1 by iV 7 (i5 7 ) ~ E 1 (dN^/dE^), then the emission rate per volume (cm -3 s~ 4 ) at cosmological time t is 

^(JSy, *) ~ "PBH W £ 7 ^( M (*)> ^ . (5-1) 

where the t dependence of M just reflects the shrinkage due to evaporation. Since the photon energy and density are 
redshifted by factors (1 + z) _1 and (1 + z)~ 3 , respectively, the present number density of photons with energy E-yg is 



<.min(t ,T) i 

n 70 (£ 7 o)= / dt(l + z)- 3 -^((l + z)£ 70 ,t) 

min(to,T) iTy 

/^•ujjI/mJ^o / dt(l + z)^(M(t),(l + z)E l0 ), 



(5.2) 



18 



where i m i n corresponds to the earliest time at which the photons freely propagate (see later) and t^pbh^o) must be 
interpreted as the current PBH number density for M > M* or the number density PBHs would have had now if they 
had not evaporated for M < M* . The photon flux is then 

c 

^=4-^70 (5-3) 

in units of s" 1 cm~ 2 sr _1 . The calculated present-day fluxes of primary and secondary photons are shown in Fig. [4j 
where the density tipbh for each M has the maximum value consistent with the observations. 




EyO [MeV] 



FIG. 4. Fluxes corresponding to the upper limit on the PBH abundance for various values of M . All PBHs produce primary 
photons but M < M* ones also produce secondary photons and this gives a stronger constraint on ft . 

Note that the highest-energy photons are associated with PBHs of mass M* . Photons from PBHs with M > M* 
are at lower energies because they are cooler, while photons from PBHs with M < M* are at lower energies because 
(although initially hotter) they are redshifted. The spectral shape differs depending on the mass M and can be 
analytically explained as follows. Holes with M > M* have a rather sharp peak which is well approximated by the 
instantaneous blackbody emission of the primary photons. On the other hand, as indicated by Eq. (|3.12l) and discussed 
in more detail in the Appendix, holes with M < M* have a high-energy n 7 o ~ E~q tail for E y o ^> ?bh/(1 + Z ( T )) 
due to their final phases of evaporation. As M increases towards M* , the primary photons start to dominate and 
generate a redshifted low- energy tail which scales as E^/q 2 jl72j |. 

Although most of the energy from both M > M* and M < M* PBHs currently appears below 100 MeV, there 
is also a higher-energy component from any PBHs with M slightly larger than M* which today have the mass m 
given by Eq. (|3.13p . The relationship between the current mass function (dn/dm) and the formation mass function 
(dN/dM) can be approximated by 



dn 

- — = min 

dm 



m \ 2 /d^\ /d7V\ 



(5.4) 



where these are interpreted as comoving number densities. The spectra are the same well above M* but one has a 
low-mass tail with dn/dm ex m 2 for m M* . In the more precise calculation presented in the Appendix, the slope 
flattens as m goes towards M* and decreases to around 1 at M w 1.5 M* . It should be stressed that the low-mass 
tail would not be present if the formation mass function were precisely monochromatic, so it is not accounted for by 
the 5 x 10 14 g curve in Fig. SJ For if all the PBHs had exactly the mass M* , they would all evaporate at exactly the 
same time. Indeed, the width of the mass function determines how much of the mass tail is present, as discussed in 
the Appendix. 
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So far we have implicitly assumed the PBHs are uniformly distributed throughout the Universe. However, if the 
PBHs evaporating at the present epoch are clustered inside galactic halos (as is most likely), then there will also 
be a Galactic background generated by PBHs with M > M* . As discussed in Sec. IVI Al the Galactic spectrum is 
essentially dominated by the solid (red) curve in Fig.|3]but with the height increased by the local density enhancement. 
While the mass tail makes a negligible contribution to the time-integrated extragalactic background, we will see that 
it is important for the Galactic background. 



Limits on P'(M) imposed by observations of the extragalactic photon background 



The origin of the diffuse x-ray and 7-ray backgrounds is thought to be distant astrophysical sources, such as blazars, 
but so far no firm consensus has bee n est ablished. The relevant observations com e fro m HEAO 1 and other balloon 
observations in the 3-500 keV range (HJ, COMPTEL in the 0.8-30 MeV range [HI], EGRET in the 30-200 MeV 
range [2251 ] . and Fermi LAT in the 200 MeV-102 GeV ra nge [226j |. In the intermediate region between HE AO 1 and 



COMPTEL, we adopt the fitting formula given in Ref. |223| . Note that the EGRET results given in Ref. [225| are 
somewhat different from the earlier ones reported in Ref. |218| . The measurements from the Fermi satellite seem more 
compatible with the earlier EGRET results. All the observations are shown in Fig. [4] 
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FIG. 5. Upper bounds on /3'(M) from the extragalactic photon background; these are conservative in the sense that no other 
contributors to the background have been subtracted. 

The upper limit on the PBH abundance is summarized in Fig. [5] and is set by the condition that the photon flux 
touches one of the upper tips of the la error bars of each observation. In what follows, we explain the form of this 
constraint analytically. We should note that the peak of the flux for each M plays an essential role in our analysis. 
We denote the energy and flux at the peak as E 7 and / , respectively, with a subscript "0" indicating current values. 

In order to analyze the spectra of photons emitted from PBHs, different treatments are needed for PBHs with initial 
masses below and above M* . We saw in Sec. IIIIBl that PBHs with M > M* can never emit secondary photons at the 
present epoch, whereas those with M < M* will do so once the peak energy of their emission becomes comparable to 
the QCD scale, which applies once M falls below M q ss 0.4M* « 2 x 10 14 g. The ratio of the peak energies of primary 
and secondary photons is inversely proportional to M because 

A ' Til 

(l + z(r))^ 1 «6T BH (xM- 1 , (l + z(r))£; c c ~^«68MeV. (5.5) 

This is consistent with the solid (red) curve in Fig. [5J which shows the ratio of E* b q / E^q derived from our numerical 
calculations. 



20 



We now use simple analytical arguments to derive the form of the primary and secondary peak fluxes. First, let 
us consider PBHs with M < Af* . Since the lifetime r for such PBHs is much shorter than t , the total number 
of primary photons emitted from a single hole is roughly M /Tbh oc Af 2 . Since the comoving number density of 
evaporating PBHs is n.pBH(T)/(l + Z ( T )) 3 oc Af~ 3 / 2 /3 from Eq. (12.41) . this being a good approximation even in the 
ACDM cosmology, the present peak flux of primary photons scales as 

/ pri oc ( Af < Af* ) . (5.6) 

(i + z(t)Y T bh 

For Af <C Af* , the temperature is always high enough to emit pions, which from Eq. (|3.10[) will produce secondary 
photons with dN^ ec / dE 7 oc Tbh oc Af -1 at the peak energy E^ oc (see Fig. [lj. Since the lifetime scales as Af 3 , this 
implies that the total number of secondary photons is proportional to Af 2 , so their peak flux scales as 

F ec oc ,., nPBH / T ^ T BH T oc M 1 ' 2 /3 (Af « AT) . (5.7) 
(1 + z(t))- 5 

This means that j sec /JP rl is constant for Af <C AT , although the energy of the photons is different. However, as the 
PBH mass increases towards AT , the temperature drops and the production of energetic (pion-producing) quarks 
and gluons becomes suppressed compared to less massive particles, so j scc /jp n decreases. For masses exceeding AT , 
pion production is exponentially suppressed and the dependence of the primary flux on M is given by 

/ pri oc npBH(to) ( Jf-) (-) oc AT~ 5 / 2 (Af > AT) . (5.8) 



V Tbh J \ t 

The dashed (green) line in Fig. [5] shows the ratio of the primary and secondary peak fluxes found by our numerical 
calculations and is well explained by the above analytic considerations. 

We now compare the predicted PBH flux with the data. The observed x-ray and 7-ray spectra roughly correspond 
to f obs oc T 7 Q 1+e ' ) , where e is small and probably lies between 0.1 (the value favored in Ref. [218j ) and 0.4 (the value 
favored in Ref. |225| ). For Af < Af* , the comparison of the primary and secondary peaks in Fig. [2] shows that the limit 
is determined by the secondary flux. The secondary peak energy is redshifted to E^ ~ m„o/[2 (1 + z(t))] oc Af 2 . 
Putting P CC (P,M) < I ohs (E* c c ) oc (T; c c )-( 1+£ ) oc Af- 2 ( 1+e ) , we can write the upper bound on f3' as 

■ M \ -5/2-2 e 



p'{M) < 3 x 10~ 27 [jfj ( M < M *) ■ (5-9) 

For Af > Af* , secondary photons are not emitted. The peak energy of the primary photons is E^q ~ Tbh oc M _1 . 
so the observations imply a limit fP ri (,3, Af) < I oha (E$) oc (T P o)~ (1+£) oc Af 1+e , which gives 



Af 



7/2+e 



/3'(M) < 4 x 10" 26 ^) (Af > AT) . (5.10) 

The precise numerical factors in the above limits have been calculated with our code but the Af dependences qual- 
itatively explain the slopes in Fig. [5] for e w 0.2-0.3. The limit bottoms out at 3 x 10~ 27 and, from Eq. (|2.5p . 
the associated limit on the density par amet er is S7pbh(AT) < 5 x 10~ 10 , which might be compared with the limit 
^pbh(AT) < 3 x 10~ 9 obtained in Ref. [222j ] . 

The limits given by Eqs. (|5.9|) and (|5.10|) are both associated with observations below 100 MeV. Although the 
PBH flux includes a high-energy T~ 2 tail for PBHs with Af < AT , this does not affect the limit on (3 because the 
tail flux falls off faster than the observed background intensity. However, we need to consider more carefully the 
constraint from observations above 100 MeV on PBHs with Af slightly larger than Af* , since these will generate the 
tail of low-mass holes described by Eq. (|5.4j) unless the initial mass function is precisely monochromatic. A related 
point is that the sudden weakening of the limit shown in Fig. [5] as Af increases through AT is entirely a consequence 
of the assumption that the PBHs have a monochromatic mass function. One would expect this discontinuity to be 
smoothed out in any realistic PBH formation scenario in which the mass function, while peaking above Af* , has a 
spread of masses such that it incorporates some Af < Af* holes. 

These nonmonochromaticity effects are discussed further in the Appendix. In particular, one must distinguish 
between a "nearly monochromatic" function centered at Af* with a finite width \i Af* , and a nearly monochromatic 
mass function centered at (1 + [i) Af* with a finite width j'Af* . In either case, one needs to decide whether the 
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definition of /3(M) , specified only for a monochromatic mass function in Eq. (|2.4[) . includes an integral over the mass 
width or finely resolves the contribution from each part of the mass band. While the second definition is more natural 
for /i 3> 1, the first one is more natural for /i <C 1 and this is what we assume here. For a mass function centered at 
M* with fi <C 1, the total number density of tail photons generated by holes of current mass m at the peak energy is 

7 tail oc UpbhCto, t Q ) - — (x m 5 /3{AQ (0 < m < M») , (5.11) 

so the condition 7 tal1 < I ohs {E^ 1 ) oc (E^ l )~^ 1+€ ^ cx m 1+e gives an upper limit on /3(M*) larger than the limit on 
P{M*) given by Eq. (|5.9|) by a factor (m/M*) c_4 . However, the limits may be comparable for /i ~ 1 because the 
mass function becomes shallower there. For a mass function centered at M > M* , the same conclusion applies; the 
limit is weaker than the constraint already shown in Fig. [5] (associated with the earlier emission of these holes below 
100 MeV) because of the m 3 effect for /i< 1 but comparable to it for fj, ~ 1. Only for v > max(l,/i) , which allows 
some PBHs to generate secondary emission, could the limit become stronger than the direct one. In this case, the 
discontinuity at M* in Fig. [5] would be smoothed out. 

Finally, we determine the mass range over which the 7-ray constraint applies. Since photons emitted at suffi- 
ciently early times cannot propagate freely, there is a minimum mass M m ; n below which the above constraint is 
inapplicable. The dominant interactions between 7 rays and the background Universe in the r elevant en ergy range 
are pair- production off hydrogen and helium nuclei and Compton scattering off electrons [1^, Il62l l227l |. The first 
dominates for E 1 > 100 MeV and is important below a redshift which is independent of E 1 ; f or th e opacity of 
0.0112 cm 2 g _1 appropriate for a 75 % hydrogen and 25 % helium mix (86|, this redshift is given by 162] 

-2/3 , q X -2/3 



l + W«H00^j y— ) , (5.12) 

with the nucleon density parameter fib and h both being normalized to modern values in the last expression. It seems 
to be entirely coincidental that this is so close to the time of recombination. The second process dominates below 
100 MeV and, in this case, z max increases with E 1 and may be somewhat smaller. The mass M m i n is then determined 
by the condition r(M m j n ) = t(z max ) . Using t(z max ) = 400 kyr and allowing for the M dependence of f(M) , we obtain 

*(^axA V3 ^ /(A^min) ^ ^ M ^ 3 x 1() 13 g . (5 13) 

We therefore extend the limit down to this mass in Fig. [5] The maximum PBH mass for which the limit is useful is 
the value at which the photon background limit crosses the density constraint Vl-ps^M) < 0.25. Using Eq. (|2.5p . this 
mass is M max ~ 7 x 10 16 g. 



VI. OTHER CONSTRAINTS FOR EVAPORATING BLACK HOLES 



The combined BBN and photon background constraints on f3'(M) are shown by the solid (red) lines in Fig. [6] In this 
section we will summarize various other constraints in the same mass range, these being indicated by the short-dashed 
(green), dotted (blue), and broken and long-dashed (red) lines in Fig. |6] Most of the green constraints have been 
studied in detail elsewhere, so they are discussed only briefly, while the blue constraints — associated with extragalactic 
antiprotons and neutrinos — are original. The nonsolid red constrain ts ar e also partly original: the long-dashed one is 
an adaptation of the Galactic 7-ray background limit found in Ref . jl07l | , while the broken one is associated with the 
damping of CMB anisotropies and is an adaptation of a previously known limit associated with decaying particles. 
The important point is that the BBN and photon background limits are the most stringent ones over almost the 
entire mass range 10 9 -10 17 g. There i s jus t a small band in the range 10 13 -10 14 g where the CMB anisotropy damping 
limit dominates. 21 cm observations [1381 ] could potentially provide a stronger constraint in some mass range around 
10 14 g, as indicated by the broken (grey) line, but such limits do not exist at present. 



A. Galactic gamma rays 



If PBHs of mass M* are clustered inside our own Galactic halo, as expected, then there should also be a Galactic 
7-ray background and, since this would be anisotropic, it should be separable from the extragalactic background. 
The ratio of the anisotropic to isotropic intensity depends on the Galactic longitude and latitude, the Galactic core 
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FIG. 6. Combined BBN and EGB limits (solid red), compared to other constraints on evaporating PBHs from LSP relics and 
CMB distortions (short-dashed green), extragalactic antiprotons and neutrinos (dotted blue), the Galactic 7-ray background, 
and CMB anisotropies (long-dashed and broken red), and the potential 21cm limit (broken grey). The density limit from the 
smallest unevaporated black holes is also shown (short-dashed green) to show where it crosses the background photon limit. 
The solid red and dotted blue limits are original to this paper. 

radius and the halo flattening. Some time ago Wright 1 106] claimed that such a halo background had been detected in 
EGRET observations between 30 MeV and 120 MeV [218j and attributed this to PBHs. His detailed fit to the data, 
subtracting various other known components, required the PBH clustering factor to be (2-12) x 10 5 hr 1 , comparable 
to that expected, and the local PBH explosion rate to be 1Z = 0.07-0.42 pc -3 yr _1 . A more recent analysis of 
EGRET data betw een 70 MeV and 150 MeV, assuming a variety of distributions for the PBHs, has been given by 
Lehoucq et al. [l07| |. In the isothermal model, which gives the most conservative limit, the Galactic 7-ray background 
requires 1Z < 0.06 pc -3 yr _1 . They claim this corresponds to 51pbh(A7*) < 2.6 x 10~ 9 , which from Eq. (I2.5|) implies 
/3'(M») < 1.4 x 10~ 26 , a factor of 5 above the extragalactic background constraint obtained in Sec. [V] Lehoucq et 
al. themselves claim that it corresponds to /3'(M„) < 1.9 x 10~ 27 but this is because they appear to use an old and 
rather inaccurate formula relating fipBH and f3 . 

It should be stressed that the Lehoucq et al. analysis does not constrain PBHs of initial mass M* because these 
no longer exist. Rather it constrains PBHs of current mass M* and, from Eq. (|3.13l) with m = M* , this corresponds 
to an initial mass of 1.26 M* . However, as shown below, this value of M does not correspond to the strongest limit 
on j3(M) . This contrasts to the situation with the extragalactic background, where the strongest constraint on f3(M) 
comes from the time-integrated contribution of the M* black holes. There would also be a Galactic contribution from 
PBHs which were slightly smaller than M* but sufficiently distant for their emitted particles to have only just reached 
us; since the light-travel time across the Galaxy is £ ga i ~ 10 5 yr, this corresponds to PBHs initially smaller than M* 
by (£ g ai/3£o) Af* ~ 10~ 5 M* , so this extra contribution is negligible. 

To examine this issue more carefully, we note that Eq. (|3.13[) implies that the emission from PBHs with initial 
mass (1 + n) M* currently peaks at an energy E w 100 (3 /i)" 1 / 3 MeV for /x < 1, which is in the range 70-150 MeV 
for 0.7 > [i > 0.08. (Note that secondary emission can be neglected in this regime, this only being important for 
/i < 0.02.) So these black holes correspond to the "tail" population discussed in Sec. IV Gl and necessarily have the 
mass function indicated by Eq. (|5.4[) . This means that the interpretation of the Galactic 7-ray limit is sensitive to 
the nonmonochromaticity in the PBH mass function, so we must distinguish between the various situations described 
in the Appendix. In particular, we need to distinguish between mass functions which center at M„ and some higher 
values of M . In the present context, we assume that the mass function is narrow, so that /9(M) can be defined as the 
integral over the entire mass width. 

The peak energy is above 150 MeV for fi < 0.08, so the 7-ray band is in the Rayleigh- Jeans region. From Eq. (|3.3[) 
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with a oc M 4 , the flux of an individual hole scales as m 3 cx /i in this regime. As clarified in the Appendix, although 
this flux must be weighted by the number density of the holes, n(m) oc m 3 oc /i for fj, -C 1, this factor is necessarily 
balanced by the ratio (AM / Am) of the mass widths at formation and now, so the limit on j3(M) scales as . 
For 0.7 > /i > 0.08, the current number flux of photons from each PBH scales as to" 1 oc ^t" 1 / 3 , so the effective 
limit on f3(M) scales as /i 1 / 3 . The observed 7-ray band enters the Wien part of the s pectr um for jj, > 0.7, so the 
limit on (3(M) weakens exponentially for M > 1.7 M* . (Although not mentioned in Rcf. [107], one could in principle 
get a stronger limit in this mass regime from observations at energies lower than 70 MeV, but we do not attempt 
this here.) Hence the largest contribution to the Galactic background and the strongest constraint on f3(M) comes 
from PBHs with M ss 1.08 M* and has the form indicated in Fig. [6] The height of the limit depends on the halo 
concentration. Although the tail holes are less numerous than the M* holes by a factor (to/M*) 3 « 3 /x , this is more 
than compensated by the halo concentration factor of 10 6 for /i > 3 x 10~ 5 . This is quite a complicated problem, so 
we have essentially adopted the model of Lehoucq et al. in determining the height. 



Since the ratio of antiprotons to protons in cosmic rays is less than 10 -4 over the energy range 100 MeV-10 GeV, 
whe reas PBH s should produce them in equal numbers, PBHs could only contribute appreciably to the antiprotons 
[88l Il08l . Il09j |. It is usually assumed that the observed antiprotons are secondary particles, produced by spallation of 
the interstellar medium by primary cosmic rays. However, the spectrum of secondary antiprotons should show a steep 
cutoff at kinetic energies below 2 GeV, whereas the spectrum of PBH antiprotons should continue down to 0.2 GeV. 
Also any primary antiproton fraction should tend to 0.5 at low energies. Both these features provide a distinctive 
signature of any PBH contribution. 

The black hole temperature must be much larger than Tbr(M^) to generate antiprotons, so the local cosmic-ray 
flux from PBHs should be dominated by the ones just entering their explosive phase at the present epoch. Such PBHs 
should be clustered inside our halo, so any charged particles emitted will have their flux enhanced relative to the 
extragalactic spectra by a factor £ which depends upon the halo concentration factor and the time for which particles 
are trapped inside the halo by the Galactic magnetic field. This time is rather uncertain and also energy dependent. 
At 100 MeV one expects roughly £ ~ 10 3 for electrons or positrons and ( ~ 10 4 for protons or antiprotons |219j . 

MacGibbon and Carr [162] originally calculated the PBH density required to explain the interstellar antiproton flux 
at 1 GeV and found a value somewhat larger than the de nsity associated with t he 7-ray limit. After the BESS balloon 
experiment measured the antiproton flux below 0.5 GeV |228f . Maki et al. |229l ] tried to fit this in the PBH scenario by 
using Monte Carlo simulations of cosmic-ray propagation. They found that the local PBH-produced antiproton flux 
was mainly due to PBHs exploding within a few kpc and used the observational data to infer a limit on the local PBH 
explosion rate of 1Z < 0.017pc -3 yr _1 . A more recent attempt to fit the antip roton data came from Barrau et al. [22j], 
who compared observations by BESS95 [22§ |. BESS98 [230| . CAPRICE [23lJ, and AMS HH with the spectrum from 
evaporating PBHs. According to their analysis, PBHs with /3'(M») w 5 x 10~ 28 would be numerous enough to explain 
the observations. This is well below the photon background limit of 3 x 10~ 27 obtained in this paper, which suggests 
that the PBHs required to explain the extragalactic photon background would overproduce antiprotons. However, 
these results are based on the assumption that the PBHs have a spherically symmetric isothermal profile with a core 
radius of 3.5 kpc. A different clustering assumption would lead to a different constraint on /3'(M*) . 

PBHs might also be detected by their antideuteron flux. Barrau et al. [2331 ] argue that AMS and GAPS would 
be able to detect the antideuterons from PBH explosions if their local density were as large as 2.6 x 10 -34 gcm~ 3 
and 1.4 x 10~ 35 gcm~ 3 , respectively. If a null result were maintained up to these levels, it would imply /3'(A/*) < 
1.5 x 10 -26 (C/10 4 )" 1 and 8.2 x 10 -28 (C/10 4 ) -1 , respectively. However, these are only potential and not actual limits. 

The most recent data come from the BESS Polar-I experiment, which flew over Antarctica in December 2004 [234| . 
The reported antiproton flux lies between that of BESS(95+97) and BESS(2000), and the p/p ratio (r sa 10" 5 ) is 
similar to that reported by BESS(95+97). Although the latter had indicated an excess flux below 400 MeV, this 
was not found in the BESS Polar-I data. However, given the magnitude of the error bars, we might still expect a 
constraint on the local PBH number density similar to that discussed above, in which case the constraint on /3(M») 
becomes 



B. Galactic antiprotons 




(6.1) 



where we have used Eq. (|2.6[) . For reasonable values of C , this is much weaker than the 7-ray background limit, but 
the value of £ is anyway very uncertain and so this limit is not shown explicitly in Fig. [S] 
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C. Extragalactic antiprotons 



Galactic antiprotons can constrain PBHs only in a very narrow range around M ~ M* since they diffusively 
propagate to the Earth on a time much shorter than the cosmic age. However, there will also be a spectrum of 
cosmic-ray antiprotons from PBHs sufficiently light that they evaporated well before the epoch of galaxy formation 
(~ 1 Gyr) . Although such pregalactic PBHs would not be clustered, the antiprotons they emitted could still be around 
today and may occasionally enter the Galaxy and be detectable at Earth, so we require that this does not happen 
frequently enough to exceed the observational limits. One might assume that the maximum mass relevant for this 
constraint corresponds to the PBHs evaporating at galaxy formation: 

1/3 / rfiijW 1/3 



M<^)' (ffl2) M.*2xl0"g. (6.2) 

This mass scale is not necessarily relevant because even the antiprotons generated inside galaxies may escape on a 
cosmological time scale and become part of the extragalactic background. However, the leakage time depends on the 
size of the cosmic-ray region, which is rather uncertain, and could well be comparable to t . We therefore adopt the 
above upper mass limit. The minimum mass relevant for this constraint is determined as follows. The mean rate of 
pp annihilations in the cosmological background, where the target p's are in hydrogen and helium nuclei, is 

rgp(t) = (<T Pp Vp) n p (t) « 2 x 1(T 22 (~j s" 1 , (6.3) 

where we have used (a pp v p ) /c s» 4 x 10 -26 cm 2 and n p (to) ~ 2 x 10~ 7 cm~ 3 . The condition that an antiproton emitted 
by a PBH with lifetime r survives until now is then 

r^-(i)(^)-^<.. 

Using Eq. Q3.7[) , we infer that the extragalactic p limit applies for 

^(■^^(/K^H^^, (6.5) 

We have seen there is no primary emission of antiprotons and the analog of Eq. (|3.10l) implies that the current number 
density from secondary emission of PBHs with mass M is 

— ""iri^w (n) {Md <") 

<i ■ ID" ;'(!/; l^LL) :■--)' cm ! . 



Since the antiproton-to-proton ratio r is required to be less than 10 5 , one finds the upper bound 

P'(M) < 3 x 10- 24 (^) (JL^ " (4 x 10 13 g < M < 2 x 10 14 g) . (6.7) 

This constraint is shown by the dotted (blue) line in Fig. [6] Although it is much weaker than the Galactic antiproton 
constraint, it is independent of the rather uncertain parameter £ . 



D. PBH explosions 

The extragalactic 7-ray background limit implies that the PBH explosion rate 1Z could be at most 10 -6 pc -3 yr _1 if 
the PBHs are uniformly distributed or 10 pc -3 yr" 1 if they are clustered inside ga lacti c halos [III [13] ■ The latter figure 
might be compared to the Lehoucq et al. Galactic 7-ray limit of 0.06pc _3 yr _1 107] and the Maki et al. antiproton 
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limit of 0.02pc~ 3 yr _1 [229j |. We now compare these limits to the direct observational constraints on the explosion 
rate. 

In the standard model of particle physics, where the number of elementary particle species never exceeds around 100, 
it has been appr eciated for a long time that the likelihood of detecting the final explosive phase of PBH evaporations 
is very low [2351 ] . However, the physics of the QCD phase transition is still uncertain and the prospects of detecting 
explo sions would be improved in less conventional particle physics models. For example, in a Hagedorn-type picture 
[2361 ], where the number o f par ticle species exponentiates at the quark-hadron temperature, the upper limit on 1Z is 
reduced to 0.05pc _3 yr _1 [2371 ] . which is comparable to the antiproton limit. 

Even without the Hagedorn effect, something dram atic is likely to occur at the QCD temperature since the number 
of species being emitted increases dramatically |l!4| . For this reason, Cline and colleagues have long argued that 
the formation of a fireball at th e QCD temper ature coul d ex plain some of the short-period 7-ray bursts (i.e. those 
with duration less than 100 ms) jTlfl, SH [239|. In Ref. [HJ they claim to find 42 BATSE candidates of this kind 
and the fact that their distribution matches the spiral arms suggests that they are Galactic. In Ref. (24lj j they claim 
that there is a class of short-period KONUS bursts which has a much harder spectrum than usua l and identify these 
with exploding PBHs. More recently they have found a further 8 candidates in the Swift data [242]. Overall they 
claim that the BATSE, KONUS, and Sw ift da ta correspond to a 4.5cr effect and that several events exhibit the time 
structure expected of PBH evaporations [243j ] . One distinctive feature of 7-ray bursts generated by PBH explosions 
woul d be a temporal delay between the high and low-energy pulses, an effect which might be detectable by Fermi 
LAT [244] . 

It is clearly important to understand how likely PBHs are to resemble 7-ray bursts fr om a theo retical perspective. 
We have seen that evaporating black holes form QCD photospheres according to Heckler [I77lll78j . A rather different 
way of producing a 7-ray burst is to assume that the outg oing charged particl es fo rm a plasma due to turbulent 
magnetic field effects at sufficiently high temperatures [l88| . MacGibbon et al. [l73| have questioned both of these 
claims and argued that the observational signatures of a cosmic or Galactic halo background of PBHs or an individual 
high-temperature black hole remain essentially those of the standard Hawking model, with little change to the detection 
probability. However, they accept that a photosphere might form for a short period around the QCD temperature, so 
perhaps the best strategy is to accept that our understanding of such effects is incomplete and focus on the empirical 
aspects of the 7-ray burst observations. 

At much higher energies, several groups have looked for 1-100 TeV photons from PBH explosions using c osmi c-ray 
detectors. However, in this case, the constraints are also strongly dependent on the theoretical model [2451 . In 
the standard model the upper limits on the explosion rate are 5 x 10 8 pc~ 3 yr _1 from the CYGNUS a rray [246| . 
8 x 10 6 pc _3 yr _1 from the Tibet array [2471 . 1 x 10 6 pc _3 yr _1 from the Whipple Cerenkov telescope 248], and 
8 x 10 8 pc~ 3 yr _1 from the Andyrchy array [249]. These limits correct for the effects of burst duration and array 
"dead time." Such limits are far weaker than the ones associated with observations at 100 MeV. They would be even 
weak er in the 45 GeV QED photosphere model advocated by Heckler, since there are then far fewer TeV particles 
250]. For example, Buga ev et al. [l89j have used Andyrchy data to obtain co nstra ints of 1 x 10 9 pc 3 yr 1 in the 
Daghigh-Kapusta model [l80Nl83 [ and 5 x 10 9 pc _3 yr _1 in the Heckler model [1771 ]. Because of the uncertainties, 
we do not show any of these limits in Fig. [Bl 



E. CMB distortion and anisotropy 

The effects of PBH evaporations on the CMB were first analyzed by Zel'dovich et al. [l92j . They pointed out 
that photons emitted sufficiently early will be completely thermalized and merely contribute to the photon-to-baryon 
ratio. The requirement that this does not exceed the observed ratio of around 10 9 leads to a limit 

^ M)<1 ° 9 (S ^ 10 "(i^i)" ( M < lo9 s)> ( 6 - 8 ) 

so only PBHs below 10 4 g could generate all of the CMB. They also noted that photons from PBHs in the range 
10 11 g < M < 10 13 g, although partially thermalized, will produce noticeable distortions in the CMB spectrum unless 

pi{M) < G£) 1 * io " 6 (i^i) 1 (i ° n g<M< 1013 g) ' (6 - 9) 

this corresponding to the condition a(M) < 1. In the intermediate mass range, there is a transition from limit (|6 .8[) 
to the much stronger limit (|6.9[) . Subsequently the form of these distortions has been analyzed in greater detail. 
If an appreciable number of photons are emitted after the freeze-out of double-Compton scattering (t > 7 x 10 6 s), 
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corresponding to M > 10 n g, the distribution of the CMB photons develops a nonzero chemical potential, leading to 
a /i distortion. On the other hand, if the photons are emitted after the freeze-out of the single-Compton scattering 
(i>3x 10 9 s), corresponding to M > 10 12 g, the distribut ion i s modified by a y distortion. These constraints were 
first calculated in t he context of decaying particle models |25l| . In the PBH context, recent calculations of Tashiro 
and Sugiyama [252] show that the CMB distortion constraints are of order f)'(M) < 1CP 21 for some range of M . The 
precise form of the constraints is shown in Fig. |6] We note that they are weaker than the BBN constraint but stronger 
than the constraint given by Eq. (|6.9[) . There are also CMB distortion constraints associated with the accretion of 
larger nonevaporating black holes and these are discussed in Sec. IVIIl 

Another constraint on PBHs evaporating after the time of recombination is associated with the damping of small- 
scale CMB anisotropics. The l imit can be obtained by modifying an equivalent calculation for decaying particles, as 
described by Zhang et ah [253j |. Their constraint can be written in the form 

log 10 C < -10.8 - 0.50 x + 0.085 x 2 + 0.0045a; 3 , x = log 10 ( 10 -i3 s -i ) > ( 6 - 10 ) 

where T is the decay rate, which corresponds to t(M) _1 in our case, and £ is equivalent to the fraction of the CDM in 
PBHs, which is simply related to (3(M) , times the fraction of the emitted energy which goes into heating the matter. 
The last factor, which includes t he e f fects of the electrons and positrons as well as the photons, will be denoted by 
/h and depends on the redshift (1621 l227j . Most of the heating will be associated with the electrons and positrons; 
they are initially degraded by inverse Compton scattering off the CMB photons but after scattering have an energy 

7 2 £cmb« 300 (^) 2 (l + z)eV«5x 10 5 (1 + *) eV w 20 (1 + z) 2 eV, (6.11) 

where 7 is the Lorentz factor and in the last expression we assume that the mass of a PBH evaporating at redshift 
z in the matter-dominated era is M ~ M* (1 + z) -1 / 2 . Since this energy is always above the ionization threshold 
for hydrogen (13.6 eV), we can assume that the heating of the electrons and positrons is efficient before reionization. 
Using Eq. (I2.5|) for /3(M) and Eq. (I3.6|) for r(M) , one can now express Eq. (|6.10|) as a limit on (3(M) . In the mass 
range of interest, the rather complicated cubic expression in M can be fitted by the approximation 

f3' {M) < 3 x IO- 30 (j^J (j^) (2-5 x 10 13 g<M<2.4x 10 14 g), (6.12) 

where /h ~ 0.1 is the fraction of emission which comes out in electrons and positrons. Here the lower mass limi t 
corresponds to black holes evaporating at recombination and the upper one to those evaporating at a redshift 6 [254l | , 
after which the ionization ensures the opacity is too low for the emitted electrons and positrons to heat the matter 
much. Equation f)6. 12[) is stronger than all the other available limits in this mass range but does not seem to have 
been pointed out before. 

Finally we note that PBH evaporations cou ld also be potentially constrained by their effect on the fo rm of the 
recombination lines in the CMB spectrum [255| . just as in the annihilating dark matter scenario [256L 1257^ . although 
we do not discuss this further here. 



F. Relic neutrinos 



Neutrinos may either be emitted directly as blackbody radiation (primaries) or they may result from the decay 
of emitted pions, leptons, neutrons, and antineutrons (secondaries). As a result, their emission spectra are similar 
to those for photons up to a normalization factor. The neutrino background can in principle constrain PBHs whose 
lifetime is more than the time of neutrino decoupling (t > Is). This corresponds to a minimum mass of M m i njt/ s» 10 9 g. 
However, the low-energy neutrinos which we have to use are poorly limited by observations. In Super-Kamiokande 
(SK), the null detection of relic v^s implies $ Pe < 1.2cm" 2 s _1 above the threshold -Ep e o = 19.3MeV [258]. As seen 
in Fig. [3 this energy corresponds to the high-energy tail for light PBHs. The constraint associated with the relic 
neutrinos is shown by the dotted (blue) line in Fig. [5] and is much weaker than the BBN and photon b ackground 
limits. Although this is an original limit, similar calculations have been made by Bugaev and collaborators [ll3lll2l| . 
However, they assume that the PBHs have a continuous mass function and they link their model with a particular 
inflationary scenario. They also consider the possibility that positrons from the s ame PBHs could explain the 511 keV 
annihilation line from center of the Galaxy, as first suggested by Okele and Rees [ll2t ] . We do not show the associated 
constraint in Fig. [5] because it is very model dependent. 
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FIG. 7. An illustration of the maximum Z/ e flux allowed by the SK limit for three different masses. For M < M* , the lower 
(upper) curves represent primary (secondary) components and the thick curves denote their sum. M > M* holes emit primary 
neutrinos only. 



G. Other relics 

Evaporating PBH should produce any other particles predicted in theories beyo nd t he standard model. The number 
ofJPBHs is therefore limited by both the abundance of stable massive particles [123] or the decay of long-lived ones 
122]. In supersymmetry (SUSY) or supergravity (SUGRA), the lightest supersymmetric particle (LSP) is stable and 
becomes a candidate for the dark matter. If LSPs are produced by the e vapo ration of the PBHs, in order not to 
exceed the observed CDM density at present, one obtains the upper bound [124J]: 

P'(M) < 10- 18 (^) (jT^gv) ~* ( M < W 11 ("W100 GeV)- 1 g) . (6.13) 

This constr aint is sho wn in Fig. [5] but depends on the mass of the LSP and is therefore subject to considerable 
uncertainty |259l l26C)j ] . In addition, unstable particles such as the gravitino or n eutralino might be produced by 
evaporatin g PB Hs. The decay of these particles into lighter ones also affects BBN [164 1261] and this gives another 
constraint [l22j : 

m ><5xlO-«(^)~ ,/2 (^)yL)- 1 (M<10» E) , (6.14) 

where Y is the limit on the number density to entropy density ratio and x§ is the fraction of the luminosity going 
into quasistable massive particles, both being normalized to reasonable values. This limit is not shown explicitly in 
Fig. |6] but it has a similar form to Eq. ()6 . 1 3|) . 



H. Reionization and 21cm signature 

The 5-year WMAP results give the optical depth as r ~ 0.1 for CMB photons emitted from the last scattering 
surface. On the other hand, recent observations of the Gunn-Peterson troughs and a 7-ray burst around z ~ 6 imply 
that reionization of the Universe occurred at z ~ 6 [254j ■ Thus PBHs cannot be so numerous that they lead to 
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reionization earlier than z ~ 6. In principle, this leads to a constraint on PBHs with 

1/3 / f(-K,\ \ 1/3 



M >M*(*-^) {^P^j «2xl0 13 g. (6.15) 



An increase in the ionization of the intergalactic medium would also produce a 21 cm signature. Mack and Wesley [136 
have shown that future observations of 21 cm radiation from high redshift neutral hydrogen could place important 
constraints on PBHs in the mass range 5 x 10 13 g < M < 10 17 g. This is essentially due to the coincidence that 
photons emitted from PBHs during 30 < z < 300 peak in the energy range in which the intergalactic medium has 
low optical depth. Any process which heats the intergalactic medium in this period will produce a signal but the 
ionizing flux of photons and electrons and positrons from PBHs would generate a distinctive feature in the 21 cm 
brightness temperature. PBHs with 5 x 10 13 g < M < 10 14 g evaporate in 30 < z < 90 and would raise the 21cm 
brightness temperature, thereby reducing the absorption seen against the CMB. PBHs with M ~ 10 14 g would raise 
the spin temperature above the CMB, so that the 21 cm appears in emission rather than absorption. PBHs with 
10 14 g < M < 10 17 g would have a less pronounced effect. The latter limit is shown in Fig. 8 of their paper and can 
be expressed in the form 

/ AT \ 7 /2 

P'(M) < 3 x 10- 29 (_j (M > 10 14 g) . (6.16) 

It bottoms out at a mass of around 10 14 g and is well below the photon background limit. The associated limits are 
shown in Fig. [fo] but only as a broken curve since they are potential rather than actual. 



VII. CONSTRAINTS ON NONEVAPORATING PBHS 



For completeness, we now review the various constraints associated with PBHs which are too large to have evapo- 
rated by now. We also include a discussion of Planck-mass relics, although these are not large — indeed they are the 
smallest conceivable objects in nature. Many of the limits assume that PBHs cluster in the Galactic halo in the same 
way as other CDM particles. In this case, Eq. (I2.5[) implies that the fraction of the halo in PBHs is related to f3'(M) 
by 

/ = 4.8 tt PBH = 4.11 xl0V(Af) (§-Y , (7.1) 

"CDM \ M Q J 

where we assume ficDM = ^-21 and this / is to be distinguished from the one appearing in the previous sections. 
Our limits on f(M) are summarized in Fig. [5] Many of them have been described elsewhere, so we only discuss these 
briefly, although we update them where appropriate. Further details ca n be found in the papers by Josan et al. jl4dl | 
(see their Table 1), Mack et al. [9(| (see their Fig. 4), and Ricotti et al. [105] (see their Fig. 9). Note that some of the 
limits are extended into the / > 1 domain, although this is obviously excluded by the independent density constraint. 
However, it is still useful to see where the limits are located since they could become stronger in the future as the 
observational data improve. 



A. Lensing constraints 



Microlensing observations of stars in th e large and small magellanic clouds probe the fraction of the Galactic halo 
in MACHOs of a certain mass range |262j . We assume that PBHs cluster in the same way as other CDM particles, so 
that Eq. (|7.1[) applies. The optical depth of the halo, tl , is defined as the probability that any given star is amplified 
by at least 1.34 at a given time. Although the relation between tl and / depends on the halo model, the so-called S 
model [89j], which is often adopted in the analysis, gives T ^ LMC ) — 4 7 x 10~ 7 / and T £ SMC ) = IAt^ MC ^ . Although 
the initial motivation for microlensing surveys was to search for brown dwarfs with 0.02 M Q < M < .08 Mq , the 
poss ibility that the halo is dominated by these objects was soon ruled out by both the MACHO 263J and EROS 
[2641 ] surveys. Instead MACHO observed 17 events apparently induced by compact objects with M ~ 0.5 M Q and 
contributing about 20% of the halo mass Recently similar claims have been mad e by the POINT- AGAPE 

collaboration, which detected 6 microlensing events in a survey of the Andromeda galaxy [265J. 

This raises the possibility that some of the halo dark matter could be PBHs formed at the QCD phase tran sition 
[ill - lisl ]. However, recent results suggest that even a 20% halo contribution of M ~ 0.5 Mq PBHs is excluded [266j |. 
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FIG. 8. Constraints on f(M) for a variety of dynamical (green), lensing (red), and astrophysical (blue) effects associated with 
large PBHs, the dominant limit for each type of effect being shown as a solid line. The effects are femtolensing and picolensing 
of gamma-ray bursts (GRB), microlensing of stars (MACHO) and quasars (QSO), millilensing of compact radio sources (RS), 
wide binary disruption (WB), globular cluster disruption (GC), dynamical friction (DF), disk heating (DH), generation of 
large-scale structure through Poisson fluctuations (LSS), accretion effects on the CMB (FIRAS, WMAP3), and gravitational 
waves (GW). The lines on the right correspond to having one PBH in the relevant volume (Galactic halo, Lyman-a forest, 
CMB anisotropy scale), the limits not applying below these lines. 



For example, the EROS-2 collaboration monitored brighter stars in a wider solid angle and obtained fewer events 
than expected. They thus obtained more s tring ent constraints on / and argued that some of the MACHO events 
were due to self-lensing or halo dumpiness [264l ] Specifically, they e xclud ed 6 x 10 -8 Mq < M < 15 M MACHOs 
from dominating the halo; combining the MACHO and EROS results |267l | extends the upper bound to 30 M . Even 
stronger constraints are obtained for an intermediate mass range. The constraints are summarized as follows: 

{1 (6 x 1(T 8 M < M < 30 M Q ) , 
0.1 (10- 6 M Q < M < 1 Mq) , (7.2) 
0.04 (10~ 3 Mq < M < 0.1 Mq) , 

The latest OGLE-II data 268, 269] yield a constraint which is consistent with this but somewhat weaker. 

Other lensing constraints come from femtolensing and picolensing of gamma-ray bursts [270j |. assuming the bursts 
are at a redshift z ~ 1, which gives / < 1 for 10~ 16 M < M < 10~ 13 M and / < 4 for 10" 12 M < M < 10" 9 M , 
respectively, microlensing of qu asars, which gives / < 1 for 10 -3 M < M < 60 M |27l| . and millilensing of compact 
radio sources, which gives [272| 

{(M/2 x 10 5 Mq)- 2 (M < 10 6 Mq) , 
0.06 (10 6 Mq <M< 10 8 Mq) , (7.3) 

(M/4 X 10 8 Mq) 2 (M > 10 8 Mq) . 

These constraints have bee n su mmarized by Josan et al. 140j and further relevant references can be found in their 
paper. Note that Hawkins [ 273j once claimed evidence for a critical density of jupiter-mass objects from observations 
of quasar microlensing and associated these with PBHs formed at the quark-hadron transition. However, the status 
of his observations is no longer clear, so this is not included in Fig. [5] 
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B. Dynamical constraints 



Biliary star sys tems with wide separation are vulnerable to disruption from encounters with PBHs or any other 
type of MACHO [274 [275]. Observations of wide binaries in the Galaxy there fore constrain the abundance of halo 
PBHs. By comparing the result of simulations with observations, Yoo et al. [2761 essentially ruled out MACHOs 
with M > 43 Mq from providing most of the halo mass. However, Quinn et al. [277} have recently performed a more 
careful analysis of the radial velocities of these binaries and found that one of the widest-separation ones is spurious. 
The resulting constraint now reads 

f(M) < I (M/4 °° M0)_1 (4 °° M ° < M ~ 103 Mq) ' (7 4) 

T{ ' [0.4 (1O 3 M <M< 1O 8 M ). { ' } 

Other dynamical constraints have been studied by Carr and Sakellariadou [278J. An argument similar to the binary 
disruption one shows that the survival of globular clusters against tidal disruption by passing PBHs gives a limit 

f(M) < { (M/3 X M0)_1 (3 X Mq < M < Mq) ' (7 5) 

[0.03 (10 6 M Q < M < 6 X 10 9 Af Q ) , 

although this depends sensitively on the mass and radius of the cluster. The limit cuts off above 6 x 10 9 Af Q because 
the encounter becomes nonimpulsive. The upper li mit o f 3 x 10 4 Mq on the mass of objects dominating the halo is 
consistent with the numerical calculations of Moore |279j | . 

Other dynamical limits come into play at higher mass scales. Halo objects will overheat the stars in the Galactic 
disc unless one has 

The upper limit of 3 x 10 6 Mq agrees with earlier calculations by Lacey and Ostriker [280], although they argued that 
black holes with 2 x 10 6 Mq could explain some fe ature s of disc heating. Positive evidence from a qualitatively similar 
effect may come from the recent claim of Totani [28 1| that elliptical galaxies are puffed up by dark halo objects of 
10 5 Mq . Another limit arises because halo objects will be dragged into the nucleus of the Galaxy by the dynamical 
friction of the spheroid stars and halo objects themselves (if they have an extended mass function), this leading to 



excessive nuclear mass unless [278 ] 



f (M/2 x 10 4 Mq)- 10 / 7 (r c /2 kpc) 2 (M < 5 x 10 5 Mq) , 
f(M) < I (Af/4 x 10 4 Mq)- 2 (r c /2 kpc) 2 (5 x 10 5 M < M < 2 x 10 6 (r c /2 kpc) Mq) , (7.7) 
[ (Af/0.1 Mq)' 1 ' 2 (M > 2 x 10 6 (r c /2 kpc) Mq) . 

The different mass regimes correspond to the drag being dominated by spheroid stars (low M), by halo objects (high 
M) and by some combination of the two (intermediate M). This limit is sensitive to the halo core radius r c but the 
dominant constituent of the halo must be smaller than 2 x 10 4 (r c /2 kpc) 1 - 4 Mq . However, there is a caveat here in 
that holes drifting into the nucleus might be ejected by the slingshot mechanism if there is already a binary black 
hole there [282| . Note that all these dynamical limits are only meaningful if there is at least one PBH per Galactic 
halo, corresponding to 

f(M) > (^—) , M halo « 3 x 10 12 Mq , (7.8) 



Mi 



lalo 



so none of them applies below this line. This is what Carr and Sakellariadou term the "incredulity limit" [278]. The 
excluded region is therefore the upper right triangle in Fig. [5] 

The effects of PBH collisions on astronomi cal o bjects — including the Earth [283| — have been a subject of long- 
standing interest [278l |. For e xample, Zhilyaev [284| has suggested that collisions with stars could produce 7-ray bursts 
and Khriplovich et al . [2851 ] have examined whether terrestrial collisions could be detected acoustically. Recently, 
Roncadelli et al. [286l j have suggested that halo PBHs could be captured and swallowed by stars in the Galactic disk. 
The stars would eventually be accreted by the holes, producing a lot of radiation and a population of subsolar black 
holes which could not be of more conventional (nonprimordial) origin. They infer a constraint / < (M/3 x 10 26 g) , 
which would give a lower limit on the mass of any objects dominating the halo. However, a careful analysis of the 
collisions of PBHs with main-sequence stars, red giant cores, white dwarfs, and neutron stars by Abramowicz et al. 
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[287( | suggests that these effects are never important. Collisions are either too rare (for M > 10 20 g) or produce too 
little power to be detected (for M < 10 20 g). Although captures of the kind envisaged by Roncadelli et al. could 
occur for M > 10 28 g, this range is already excluded by the microlensing constraints. We therefore do not show any 
collisional constraints in Fig. [5J 

Each of these dynamical constraints is subject to certain provisos but it is interesting that they all correspond to an 
upper limit on the mass of the objects which dominate the halo in the range 500-2 x 10 4 M & , the binary disruption 
limit being the strongest. This is particularly relevant for constraining models, such as that proposed by Frampton 
et al. [40J, in which the dark matter is postulated to c omp r ise P BHs in this mass range. However, apart from the 
Galactic disc and elliptical heating arguments of Refs. [280l l28lj . it must be stressed that none of these dynamical 
effects gives positive evidence for MACHOs. Furthermo re, n one of them requires the MACHOs to be PBHs. Indeed, 



they could equall y w ell be clusters of sm aller objects [2881 ]; this is particularly pertinent in light of the claim by 



Dokuchaev et al. [71| and Chisholm [289] that PBHs could form in tight clusters, giving an overdensity well in excess 
of that provided by the halo concentration alone. 

C. Large-scale structure constraints 

Sufficiently large PBHs could have important consequences for large-scale structure formation because of the Poisson 
fluctuations in their number density. This was first pointed out by Meszaros [9 i f. alt houg h he overestimated the 
magnitude of the effect (92j. Subsequently, various authors have invoked this effect [90l I93l495l ] In particular, Afshordi 
et al. 95] use observations of the Lyman-a forest to obtain an upper limit of about 10 4 M Q on the mass of any PBHs 
which provide the CDM. Although this conclusion is based on detailed simulations, one can understand their result 
heuristically and extend it to the case where the PBHs only provide a fraction f(M) of the dark matter by noting 
that the initial density perturbation associated with the Poisson fluctuations in the number of PBHs on a mass scale 
M Lya is 

Since these fluctuations can grow between the redshift of CDM domination (z oq ~ 4000) and the redshift at which 
Lyman-a clouds are observed (zLya ~ 4) by a factor z cq / z\^ ya ~ 10 3 , imposing the condition that S(M) has not grown 
to more than order unity (as required for Lyman-a clouds) leads to the constraint shown in Fig. [8j 

f{M)< (wk)'^ (7 - 10) 

This form of the statistical clustering limit was also derived in Ref. [92j and is shown in Fig. 5 of Ref. [9(| ■ Note that 
this limit requires at least one PBH per Lyman-a mass, which corresponds to 

f(M) > (j^J > ^Lya « 10 10 M . (7.11) 



The data from SDSS are now more extensive [2901 ] . thereby reducing the PBH contributi on to the fluctuations (i.e. 
reducing the value of 6 at z = ^Lya)- Indeed, the limiting mass might be reduced to 10 3 M [29l| . which is comparable 
to the strongest dynamical limits discussed above, although we do not use this in Fig. [8] A similar effect can allow 
clusters of large PBHs to evolve into the supermassive black holes thought to reside in galactic nuclei [96l - [98| ; if one 
replaces M-^ yoi with 10 s M and z\^ ya with 10 in the above analysis, one reduces the limiting mass in Eq. (|7.10[) to 
600 Mq . However, this limit is not shown in Fig. [5] 

D. Accretion constraints 

There are good reasons for believing that PBHs cannot grow very much during the radiation-dominated era. 
Although a simple Bondi-type argument suggests that they could grow as fast as the horizon Q, this does not 
account for the backgro und cosmological expansion and a fully relativistic calculation shows that such self-similar 
growth is impossible [3, |292| . |293| ] . Consequently there is very little growth at all. The only exception might be if 
the Universe were dominated by a "dark energy" fluid with p < — (1/3) pc 2 , as in the quintessence scenario, since 
self-similar black hole solutions do exist in this situation [2911296]. This may support the claim of Bean and Magueijo 
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[99J that intermediate mass PBHs might accrete quintessence efficiently enough to evolve into the supermassive black 
holes in galactic nuclei. 

Even if PBHs cannot accrete appreciably in the radiation-dominated era before decoupling, massive ones might still 
do so in the period after decoupling and the Bondi-type analysis should then apply. The associated accretion and 
emis sion of radiation could have a profound effect on the thermal history of the U niver se, as first analyzed by Carr 
[297| . Recently this possibility has been investigated in great detail by Ricotti et al. j 1051 ] , who study the effect of such 
PBHs on the ionization and temperature evolution of the Universe. They point out that the x rays emitted by gas 
accretion onto the black holes could produce measurable effects in the spectrum and anisotropies of the CMB. Using 
FIRAS data to constrain the first and WMAP data to constrain the second, they impro ve th e constraints on f(M) 
by several orders of magnitude for M > 0.1 M . These limits are shown in Fig. 9 of Rcf. 105], although their precise 
form clearly depends on details of the accretion model. They should also be updated in principle to comply with the 
WMAP 7-year results. Note that this limit requires at least one PBH on the scale associated with the relevant CMB 
anisotropies; for I = 100 modes, this corresponds to 

f(M) > ( -J*— ) , M l=lm « 10 16 M . (7.12) 

In related work, Mack et al. [9(| have considered the growth of large PBHs through the capture of dark matter halos 
and suggested that their accretion could give rise to ultraluminous x-ray sources. The latter possibility has also been 
explored by Kawaguchi et al. [Hj|. 



E. Gravitational-wave constraints 

A population of massive PBHs would be expected to generate a background of gravitational waves [l03. This 
would be especially interesting i f there w ere a population of binary halo PBHs coalescing at the present epoch due 
to gravitatio nal r adiation losses jlOll . Il02| . If such waves were detected, it would pro vide a unique probe of the halo 
distribution [l03j | . The LIGO data already place weak constraints on such scenarios 298] . Howev er, a diffe rent type 
of gravitational- wave constraint on f3(M) has recently been pointed out by Saito and Yokoyama [299l |300j. This is 
because the second-order tensor perturbations generated by the scalar perturbations which produce the PBHs are 
surprisingly large. The associated frequency is around 10 -8 (M/10 3 M Q ) Hz, while the limit on j3(M) just relates to 
the amplitude of the density fluctuations at the horizon epoch and is formally of order 10 -52 . The limit from pulsar 
timing data is shown in Fig. [8] and already excludes PBHs with 10 2 Mq < M < 10 4 M© from providing an appreciable 
amount of dark matter. The potential limits from LIGO, VIRGO, and BBO are not shown exp licitly but cover the 
mass range down to 10 20 g. This effect has also been noted by several other authors |30ll . l302j [. These limits apply 
if and only if the PBHs are generated by super-Hubble scale fluctuations, like the ones created during inflation. We 
also note that the limiting value of j3 depends on the fluctuations being Gaussian. Although this is questiona ble in 
the context of the large- amplitude fluctuations relevant to PBH formation, the studies in Refs. [25|, 1 144 Il45| show 
that non-Gaussian effects are not expected to be large. Gravitational-wave detectors in space, like LISA, might 
also constrain PBHs indirectly. This is because LISA could detect isola ted PBH s in the mass range 10 14 -10 20 g by 
measuring the gravitational impulse induced by any nearby passing one [303l |304| . 



F. Planck-mass relic constraints 



If PBH evaporations leave sta ble Planck-mass relics, these might also contribute to the dark matter. This was 



first pointed out b y MacGib bon [128l ] and has subsequently been explored in the context of inflationary scenarios by 
several authors (8l. Il29t — 135] . If the relics have a mass nMp\ , where Mpi is the Planck mass, and reheating occurs at 
a temperature Tr , then the requirement that they have less than the critical density implies Q 



for the mass range 



M \ 3/2 

0'(M) < 8 x 10- 28 k- 1 [ — ) (7.13) 



2 < — < 10 11 k 2 '" . (7.14) 
TpJ Mpi V ' 



Note that we would now require the density to be less than ficDM ~ 0.25, which strengthens the limit by a factor of 
4. The lower mass limit arises because PBHs generated before reheating are diluted exponentially. The upper mass 
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limit arises because PBHs larger than this dominate the total density before they evaporate, in which case the final 
cosmological photon-to-b aryo n ratio is determined by the baryon asymmetry associated with their emission. Recently 
Alexander and Meszaros jl36l | have advocated an extended inflationary scenario in which evaporating PBHs naturally 
generate the dark matter, the entropy, and the baryon asymmetry of the Universe. This triple coincidence applies 
providing inflation ends at t ~ 10 -23 s, so that the PBHs have an initial mass M ~ 10 6 g. This just corresponds to 
the upper limit indicated in Eq. (|7.14j) . which explains one of the coincidences. The other coincidence involves the 
baryon asymmetry generated in the evaporations. It should be stressed that the limit (|7.13l) still applies even if there 
is no inflationary period but then extends all the way down to the Planck mass. 



VIII. CONCLUSIONS 

All the limits considered in this paper are brought together in a master f3'(M) diagram in Fig. |9] In particular, 
the constraints on f(M) discussed in the previous section have been converted into limits on (i'(M) using Eq. (17.11) . 
We also include the relic limit associated with Eq. (|7.13p — with the broken line to the left applying if there is no 
inflation — and the entropy limit associated with Eq. (|6.8[) . The latter is also shown broken since it is much weaker 
than the LSP constraint, albeit more secure. Most of the limits are associated with various caveats, but where 
they are reasonably firm, only the dominant one is indicated for each value of M . Nevertheless, we include several 
overlapping ones at high masses. Figure [9] covers the entire mass range from l-10 50 g and involves a wide variety of 
physical effects. This reflects the fact that PBHs provide a unique probe of the early Universe, gravitational collapse, 
high-energy physics, and quantum gravity. In particular, they can probe scales and epochs inaccessible by any other 
type of cosmological observation. 




10 20 30 40 50 

log 10 (M/g) 

FIG. 9. Master fi'(M) constraints diagram for the mass range 1-10 50 g, the acronyms being specified in the caption to Fig. [8] 

Although none of the effects discussed in this paper provides positive evidence for PBHs, Fig. |H] illustrates that 
even the nondetection of PBHs allows one to infer important constraints on the early Universe. In particular, the 
limits on f3(M) can be used to constrain all the PBH formation mechanisms described in Sec. HI Thus, for example, 
they constrain models involving inflation, a dustlike phase, and the collapse of cosmic strings or domain walls. They 
also restrict the form of the primordial inhomogeneities (whatever their source) and their possible non-Gaussianity. 
Finally, they constrain less conventional models, such as those involving a variable gravitational constant or extra 
dimensions. However, it must be emphasized that the form of the /3(M) limits may itself change in such models, so 
it is not just a matter of applying the form of the limits derived in this paper directly. These issues are too broad to 
address here but they provide much scope for future work. 
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Appendix: Link between mass tail, energy tail, and nonmonochromaticity of PBH mass function 

The high-energy tail generated by PBH evaporations has been studied in detail by MacGibbon |l72j . so we only 
summarize the argument briefly here. We confine attention to the photon tail for simplicity. The form of this tail can 
be derived by integrating Eq. (13.21) over time (i.e. M) for fixed E : 

dE J Mi dE dM J e EM -l ' V ; 

where we have used the relation dM/dt oc M -2 , a s (E,M) is the absorption cross section appearing in Eq. (13.31) . 
which scales as M 2 for E > M~ x and E 2 M 4 for E < M" 1 (l69j . M ; is the initial mass of the hole and (for the 
current discussion) M is the evolving mass. We are here neglecting redshift effects, because most emission occurs at 
fixed redshift for a given value of M{ . However, this could be accounted for by replacing E by the present photon 
energy Eq = E (1 + Zcvap)^ 1 • We also drop the suffix j on E . 

Let us first assume Mj > M q s» 0.4M, , so that secondary emission is initially unimportant. For E < Mr 1 , the 
mass integral just involves the Rayleigh- Jeans part of the spectrum, so it is dominated by the upper limit and scales as 
Mf E , leading to dN 1 /dE oc E 3 . For E > Mr 1 , the exponential term cuts the integral off above a mass M ~ E^ 1 , 
so the integral scales as E~ 5 , leading to dN 1 /dE oc E~ 3 . The time-integrated spectrum of photons from a PBH with 
Mj > M q can therefore be expressed as 



diV-r pc \E 3 (E < M i 
dE [E- 3 (Air 1 < e < A QC d) 



(A.2) 



Once se cond ary emission becomes important, as is always the case for M; < M q , the analysis of MacGibbon and 
Webber Il6lll shows that 



dA 7 ^ (E- 1 (m^<E < M- 1 ) , 



dE 



E 2 M 2 e -BM \ E >M- 1 ) 



(A.3) 



The form around the peak at m n reflects the low-energy fragmentation function, which is roughly Gaussian, although 
we do not give this explicitly. One now uses this expression in Eq. (|A. 1|) and integrates over M . For E < Mr 1 , 
the integral is dominated by the emission when M ~ M; , so dN 1 /dE scales as E~ x . For E > Mr 1 , the integral is 
dominated by holes with M ~ E^ 1 , so dN 7 /dE scales as E~ 3 . Therefore the E~ 3 tail given by Eq. (|A.2[) simply 
extends into the E > Aqcd regime for Mj > M q . However, for M; < M q , we obtain 




(Mr 1 >E> Aqcd) , (A , 

(E > Mr 1 ) , 



where the form around is not shown explicitly but again reflects the low-energy fragmentation function. The 
important qualitative point is that one has the same high-energy E~ 3 tail as before, although there is now an 
intermediate E^ 1 regime. 

There is also a contribution to the photon background above 100 MeV from the current remnants of PBHs which 
were initially slightly larger than M* . This is distinct from the high-energy tail discussed above, although (as we will 
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see) the two contributions are related. This can be understood as a consequence of Eq. (|3.13[) , since this implies that 
the current mass function {An/ Am) is related to the formation mass function (AN /AM) by 

An _ frn_\ 2 ( 1 \ 2 / AN \ _ / mV /dA 



both mass functions being comoving. The first expression is exact, with /j,(m) being implicitly determined by 
Eq. (|3.13[) . while the second expression applies for fi <C 1- In the latter case, m ss (3 ^t) 1 ^ 3 M» and the integrated 
comoving number density of holes with mass between and m can be approximated by 

1 / m \ 3 /dA\ 1 / m x ; 



" (<m) ^3UJ M^A^UJ (m<<M * h (A - 6) 

where we have used Eq. (I2.9[) . A correction would be required if the spectrum did not extend all the way down to 
zero mass (i.e. if the initial spectrum did not extend all the way down to M„). We describe this low- mass part of the 
present spectrum as the "mass tail." Note that An/ Am « AN /AM for /i S> 1 and this reflects the fact that m « M 
once [i goes much above 1. [For example, Eq. (|3.13[) implies the ratio m/M is (26/27) 1 / 3 for fi — 2 and (63/64) 1//3 
for u — 3.] There is no simple analytic expression for n(m) in the intermediate case (/i ~ 1) but the overall current 
mass function can be approximated by 



dn ' 
dm 



m\ f AN\ f AN 
WJ \AMj J [dM 



(A.7) 



In a more precise calculation, Eqs. (|A.5p and (|3.13[) imply that the slope of the mass function decreases as m 
approaches M* , with the exponent of m being 2/(l+/i) 3 . The relationship between An /Am and AN /AM is represented 
qualitatively in Fig. [101 although it should be noted that this assumes AN/ AM is a decreasing function of M . 
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FIG. 10. This shows the relationship between the PBH mass functions at formation and currently. For a nearly monochromatic 
initial mass function centered at M» , the extent of the current low-mass tail depends on the initial mass width. It covers only 
low masses if AM <C M* [case (1A)] but the entire range from to M* if AM ~ M 4 [case (IB)]. For a nearly monochromatic 
mass function centered above M* , we have either the low end of the tail [case (2A)], part of the upper end [case (2B)] or the 
entire tail [case (2C)]. 



It should be stressed that the low-mass tail would not be present if the formation mass function were precisely 
monochromatic. For if all the PBHs had exactly the mass M* , they would all evaporate at exactly the same time. 
Indeed, the width of the mass function determines how much of the mass tail is present. In order to understand these 
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nonmonochromaticity effects, we must distinguish between two different situations: (1) a "nearly monochromatic" 
function centered at M* but with a finite width /zM*, so that it extends up to (1 + fi) M* ; and (2) a "nearly 
monochromatic" function centered at (1 + fx) M* but with a finite width v M* (where \i may be different from v). 
In the first situation, there are two subcases: for \i -C 1 [case (1A)], we only have the lower part of the m 2 tail (i.e. 
the spectrum does not extend all the way up to M*); for fj, ~ 1 [case (IB)], we have the entire tail. In the second 
situation, there are three subcases: for /x -C 1 and v <C 1 [case (2A)], we only have part of the lower mass tail; for 
/i ~ 1 and i/« 1 [case (2B)], we only have part of the tail close to M*; for v > max(l,/z) [case (2C)], we have the 
entire tail. Clearly cases (2A) and (2B) are most consistent with the rationale of the present analysis, which is to 
assume a monochromatic mass function. The different possible forms of the current mass function dn /dm in these 
subcases are illustrated in Fig. [TU] 

Since the photon production rate of an individual hole is _/V 7 oc m" 1 , the instantaneous flux from the tail population 
is / oc n(m) mT 1 cx m 2 oc E~ 2 , which is reminiscent of the high-energy tail from the PBHs with M < M* . However, 
the connection between the mass tail for M > M* and energy tail for M < M* is a subtle one, which requires some 
clarification. All PBHs generate an E~ 3 energy tail eventually but only those with M < M* produce one by the 
present epoch and PBHs with mass slightly above M* do not produce the entire energy tail because they have still not 
completed their evaporation (i.e. the highest-energy part is missing). But it is precisely these unevaporated remnants 
which provide the mass tail, so the energy and mass tails are complementary in the sense that the latter becomes 
important just as the former becomes unimportant. 

Finally, it should be stressed that the definition of f3(M) requires some care for an extended mass function, the 
definition given by Eq. (|2.4[) only applying in the monochromatic case. One needs to decide whether to integrate 
dN/dM over the entire mass band or over smaller ranges within this band. While the second definition is more 
natural for /i> 1, the first one is more natural for fi <C 1 and this is what we assume in this paper. In discussing the 
value of fi(M) for M ~ M* , one must also distinguish between the width of the mass band at formation, AM , and 
its width at the current epoch, Am . For a narrow mass function, the PBH collapse fraction at formation is 

« M - )a (^), aM =(^) am - (A - 8) 

Although dn/dm is less than dN/dM by a factor (m/M*) 2 for m <C M* , this is necessarily balanced by the factor 
(AM / Am) , since the PBH number density itself is preserved prior to complete evaporation. So there is no extra 
(m/M*) 3 factor in any constraint on (3(M) . However, the current integrated number density of black holes between 
and m is less the number density of the original M* black holes by a factor (m/M») 3 because this corresponds to 
PBHs with initial mass in some small range above M* . 
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